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Abstract 

An amplitude equation for an unstable mode in a collisionless plasma is de- 
rived from the dynamics on the unstable manifold of the equilibrium Fq(v). 
The mode eigenvalue arises from a simple zero of the dielectric ej.(z); as 
the linear growth rate 7 vanishes, the eigenvalue merges with the continu- 
ous spectrum on the imaginary axis and disappears. The evolution of the 
mode amplitude p(t) is studied using an expansion in p. As 7 — > + , the 
expansion coefficients diverge, but these singularities are absorbed by rescal- 
ing the amplitude: p(t) = 7 2 r(7i). This renders the theory finite and also 
indicates that the electric field exhibits trapping scaling E ~ 7 2 . These sin- 
gularities and scalings are independent of the specific Fq(v) considered. The 
asymptotic dynamics of r(r) can depend on Fq only through expi£ where 
dek/dz = \e' k \ exp — Similar results also hold for the electric field and 

distribution function. 
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I. INTRODUCTION 



The evolution of an unstable electrostatic mode is a fundamental problem in collision- 
less plasma theory. Although quite idealized, this evolution involves many features that are 
essential to more complicated and realistic problems, in particular there is a singular inter- 
action between the wave and resonant particles. This resonance drives the initial growth of 
the unstable linear mode, and then trapping of resonant particles by the finite amplitude 
wave marks the onset of strong nonlinear effects that saturate the instability. These nonlin- 
ear effects are difficult to treat analytically while maintaining the self-consistent relationship 
between electric field and particles, and calculations on this problem have emphasized spe- 
cial regimes which allow simplifying approximations, e.g. a "bump on tail" distribution or 
instability driven by a small cold beam. [l| [TTJ] 



In this paper, I describe a new approach which simplifies the problem by restricting 
attention to the dynamics occuring on the unstable manifold of the equilibrium. Physi- 
cally this restriction means I consider initial conditions in which only the unstable modes 
are excited, rather than allowing arbitrary initial conditions comprised of all linear modes. 
Mathematically the unstable manifold is finite-dimensional and this reduction in dimension 
provides a considerable simplification. 

The restriction on initial conditions is compensated by the freedom from inessential 
assumptions about the equilibrium F (v, ji). Here \i denotes any parameters such as density 
or temperature that determine the properties of Fq; it is not necessary to make a specific 
choice for /i, rather we let F (f,/i c ) denote the critical equilibrium. For /i < /i c , F (t>,/i) is 
linearly stable and for fi > fi c there is an unstable mode (or modes) with linear growth rate 
7. The limit fi — > fi c from the unstable regime will usually be denoted 7 — > + . Thus I am 
able to give a unified treatment of instabilities in beam-plasma systems with warm or cold 
beams as well as two-stream instabilities for equilibria with counterstreaming components 
of equal density. 

An additional motivation for this approach is the possibility that the dynamics on the 
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manifold will exhibit universal features of the instability. This hope arises from experience 
with simpler bifurcations in dissipative dynamical systems. Consider, for example, a Hopf 
bifurcation [0,0 in which an equilibrium loses stability as an isolated conjugate pair of 
eigenvalues (A, A*) cross the imaginary axis while all other modes remain stable. In this 
situation, there is a two-dimensional unstable manifold associated with the unstable modes, 
and the isolation of critical eigenvalues permits the dynamics near the equilibrium to be 
rigorously reduced to the two-dimensional dynamics on the unstable manifold. In polar 
variables (p, 0) , this reduced two-dimensional system has the form 

p = 7P + a lP 3 + a 2 p 5 + 0{p r ) (1) 
= lu + a[p 2 + a' 2 p 4 + 0(p 6 ) (2) 

where A[t) = p{t) e~ ld ^ is the amplitude of the unstable linear mode and A = 7 — iuj 
is the critical eigenvalue; the evolution of pit) decouples from the phase 6{t) so that the 
dynamics can be easily analyzed. Provided the cubic coefficient a\ is non-zero at the onset 
of instability, these equations have a universal structure in the limit of weak instability 
7 —>■ + . More specifically, by setting p(t) = y/j r( r yt) we obtain from ([[]) 

dv 

— = r + air 3 + 7 a 2 r 5 + C( 7 2 ) (3) 
ar 

where r = jt. As 7 — > + , the terms of higher order in r vanish leaving f = r + air 3 , an 
asymptotic equation of universal form reflecting the specific problem under consideration 
only through the coefficient a%. In this way, the unstable manifold dynamics for Hopf 
bifurcation reveals the slow time scale r and, more interestingly, a universal scaling behavior 
p ~ yFy for the mode amplitude. 

There are many key differences between such a simple dissipative bifurcation and the 
bifurcation arising from the appearance of an unstable mode in the linear spectrum of 
a Vlasov equilibrium. The Vlasov equation is a Hamiltonian dynamical system and the 
spectrum for a stable equilibrium is pure imaginary. |T^-|TB| The eigenvalues describing 
the unstable modes are not isolated at the onset of instability, in fact they appear in the 

3 



spectrum for the first time at onset and are embedded in the continuous spectrum on the 



imaginary axis. |I7| Furthermore, without dissipation, one does not expect the unstable 
manifold associated with the instability to be attracting. 

Despite these important differences, it is feasible to adapt the derivation of ([I]) - (0) to 
the bifurcation of an electrostatic mode and to obtain the corresponding equations for the 
dynamics on the unstable manifold. The 7 — > + limit of these equations is remarkable 
and provides a striking contrast with the familiar limiting behavior in Hopf bifurcation. [18| 



For an electrostatic mode, the coefficients (aj,a'j) are singular as 7 — > + with asymptotic 
form 7~( 4 -? -1 ). These divergences are not unphysical, rather they signal the presence of a 
different scaling from that characterizing a Hopf bifurcation; by setting p(t) = 7 2 r(7t) one 
can absorb the singular behavior at every order leaving rescaled equations for r that are well 
behaved as 7 — > + . Unlike the rescaled equation for Hopf bifurcation (|3]), in the Vlasov 
case the terms that are higher order in r are not higher order in 7, rather at each order 
bj r 2j the rescaled coefficient bj = aj is order unity as 7 — >• + . Hence the equation 

for dr/dr does not truncate, and the 7 — > + equation retains an infinite set of terms. 
The first two coefficients have been calculated and shown to be independent of F Q (v,p c ) as 
7 — > + with values 61 = —1/4 and 62 — 13/64. The coefficients at higher order are not 
explicitly determined, however one can prove that at each order there is a universal function 
Qj, independent of F (v,p c ), such that as 7 — > + 

bj = ReQjie^) (4) 

where the phase = e' k * /e' k is defined in terms of the 7 — > + limit of the derivative of the 
dielectric function e^. The identification of as a scaling variable for the 7 — >• + regime 
which captures any remaining dependence on the underlying equilibrium is a novel result of 
this work. 

The scaling p(t) ~ 7 2 for the amplitude of the unstable mode implies that the wave 
electric field follows the so-called trapping scaling E ~ 7 2 in the limit 7 — > + . The termi- 
nology arises from the equivalent scaling u>j, ~ 7 between the bounce frequency of trapped 



particles and the growth rate (u% = ekE^/m). Trapping scaling has been a characteristic 



feature of previous numerical simulations 0|| as well as recent experiments |19[; however 
the theoretical results on the scaling of the saturated electric field are divided between theo- 
ries predicting "Hopf scaling" E ~ ^/j |1|J^,|1^,|11|] , and theories predicting trapping scaling. 



>|. |20fl Perhaps the best known analysis leading to a prediction of Hopf scaling is a con- 
troversal paper by Simon and Rosenbluth. @ Their work treats a one mode bump on tail 
instability by perturbatively expanding the Vlasov equation and seeking a time-periodic 
nonlinear solution. The perturbation theory leads to singular results and the final expres- 
sions are rendered finite by prescribing a regularization procedure. Subsequent perturbation 
theories have encountered comparable difficulties and proposed similar prescriptions. UT0|JTTJ| 
The approach I develop in this work differs crucially from these investigations in the 
interpretation and treatment of the singular behavior of the expansion at 7 = 0. The 
derivation of the amplitude dynamics on the unstable manifold gives nonlinear coefficients 
as integrals over velocity, and as 7 — > + these integrals diverge due to pinching singularities 
that develop at the phase velocity of the mode. As already mentioned, these divergences can 
be absorbed by simply rescaling the wave amplitude and this rescaling reveals an electric 
field that follows the trapping scaling. If, on the other hand, one were to regularize the 
integrals by somehow discarding the divergent part then the resulting equations for the 
mode amplitude would indeed scale as in dissipative Hopf bifurcation. This, in essence, is 
the step taken in the theories that predict Hopf scaling for the electric field. 

The theory of invariant manifolds for equilibria of infinite-dimensional dynamical sys- 
tems, such as partial differential equations, has been developed extensively in recent years 
with rigorous results establishing the existence and properties of these structures for vari- 
ous classes of evolution equations. Examples utilizing a variety of techniques are found in 
pl^p6[, and there is an introductory review by Vanderbauwhede and Iooss. |27j However, 



this theory does not yet treat equations such as the Vlasov-Poisson system and it seems 
to be an open problem to rigorously construct invariant manifolds for Vlasov equilibria. In 
this paper, unstable manifolds serve a heuristic role by motivating certain procedures for 



constructing the expansions in the mode amplitude. Although these same expansions could 
certainly be set up without mentioning the manifolds, the dynamical systems viewpoint 
does seem to bring additional insight. Future development of a rigorous invariant manifold 
theory for the Vlasov equation may provide mathematical justification for these expansions. 
The present analysis is solely concerned with understanding the properties of the instability 
as represented by the amplitude expansions. 

The remainder of the Introduction is devoted to defining our notation and summarizing 
relevant well known facts concerning the spectrum and eigenf unctions of the Vlasov-Poisson 
equation. In Section II the description of the unstable manifold is briefly reviewed and the 
equations necessary to obtain the dynamics on the manifold are derived. These equations are 
solved using power series in the amplitude of the unstable mode in Section III and a detailed 
analysis of the lowest order term in this expansion reveals the divergence mentioned above. 
This divergence is shown to imply the trapping scaling for the mode amplitude. In Sections 
IV- VI, the structure of the expansions is examined to all orders. The increasing strength 
of the divergences is calculated and a detailed analysis is made of how the dynamics on the 
unstable manifold depends on the critical equilibrium F (v,fj, c ) in the 7 — > + limit. The 
mode amplitude dynamics, the electric field, and to a large extent the distribution function 
depend on F (v , u c ) only through the derivative of the dielectric function e' k . This conclusion 
indicates a degree of universality to the dynamics of a weakly unstable electrostatic mode 
that has not been previously appreciated. 



A. Notation 

For a neutral collisionless plasma with a fixed ion density no, the electron distribution 
function F(x,v,t) and the electrostatic potential &(x,t) satisfy the dimensionless Vlasov- 
Poisson equations (in one dimension) 

— + v— + — — -0 (5) 

dt dx dx dv 
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The plasma length is L and periodic boundary conditions are assumed with normalization 



/L/2 roo 
dx / dvF(x,v,t) = L. (7) 
-L/2 J-oo 



Here x, t and v are measured in units of u/u p , u p 1 and u, respectively, where u is a chosen 
velocity scale and uj^ = 47re 2 n /m. The electron charge and mass are — e and m and the 
ions are singly charged. The dimensionless distribution function and potential are measured 
in units of w -1 and mu 2 /e respectively. 

An inner product between two functions Gi(x,v) and (^(x, v) is defined by 



(G 1 ,G 2 )= f dx <G 1 ,G 2 > (8) 

J -L/2 

where < Gi, G 2 > denotes the integration over v alone: 

/oo 
dvG 1 (x,v)*G 2 (x,v). (9) 
-oo 

With periodic boundary conditions the allowed wavevectors are multiples of k c = 2tt/L, and 
the Fourier expansion of G(x, v) will be written as 

oo 

G{x,v)= E e ikk ° x G k (v). (10) 

k=— oo 

Let Fo(v,n) denote a parametrized family of equilibria, normalized to unit density 

/oo 
dvF (v,fi) = l, (11) 
-oo 

and re-express the distribution function relative to this family F(x,v,t) = F (v,fi) + 
f(x,v,t), then / satisfies 

W +v W + ?±f^ + ?±?l = (12) 
dt dx dx dv dx dv 



dx 2 



dv f(x,v,t) (13) 



and 



L/2 poo 

dx / dv f(x,v,t) = 0. (14) 

-L/2 J—oo 



With the Fourier series for /, equations flT2|) - ( |T5D can be combined 

5/ 



Cf + M(f) (15) 



where the linear operator is defined by 



oo 



Cf= E e 1 **-* (16) 

fc=— oo 



fc = 

(L h f k )(v) = { r n (17) 



vfk{v) + k~ 2 T](v, fi) J? dv' f k (v') k^O 



with 



OF 



and 



W) = £ E E 'y^ 1 /" ^7KA (19) 



oo 



oo 



c fc=— oo {=— oo 

Here and below a primed summation omits the Z = term. 

Symmetries of the model @ - @ and the equilibrium F (v,fi) are important qualita- 
tive features of the problem. Spatial translation, T a : (x,v) — > (x + a, t>), and reflection, 
k : (x, f) — > (— x, — act as operators on the distribution function in the usual way: if a 
denotes an arbitrary transformation then (a ■ f)(x,v) = /(a™ 1 • (x,v)). The operators £ 
and M commute with T a due to the spatial homogeneity of F , and if F (v, /i) = Fq(—v, /z), 
then £ and A/" also commute with the reflection operator k. Together T a and k generate the 
symmetry group of the circle 0(2) and without k the symmetry drops to SO (2). 

The spectral theory for £ is well established, and the needed results are simply recalled 
to establish the notation. flTT| , [2"8"| , P§| The eigenvalues A = —ikk c z of £ are determined by the 
roots Ak(z,fj) = of the "spectral function", 
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v — z 



(20) 



Unless it is necessary to manipulate the parameter dependence of A^(z, /i), the argument fi 
will be suppressed. If only translation symmetry is present, these eigenvalues are generically 
complex; when the equilibrium is also reflection-symmetric then real eigenvalues can arise, 
for example in a two-stream instability. [B(J 

The spectral function is analytic in the upper and lower half planes with a branch cut 
on the real axis along the support of rj. For z = r ± ie, the discontinuity across the cut is 
given by the Plemelj formula [|3TJj 



lim Afc(r ± ie) = 1 + — 



A' 2 



p.v. r 

oo v — r 



± iirr)(r, //). 



(21) 



The analytic continuation of A&(z) from Im z > to Im z < yields the dielectric function 
£k( z ) |32]| ) defined in the usual way via the Landau contour. |33[] Since our analysis focuses 
on the regime Im z > 0, the notations Ak(z) and e^(z) are interchangeable in the subsequent 
discussion. 

The branch cut of A&(z) corresponds to a continuous spectrum for L on the imaginary 
axis. As fi varies the roots of A&(z) typically vary; in particular, roots can appear or 
disappear through the branch cut. The appearance of a root at the cut corresponds to the 
birth of eigenvalues embedded in the continuous spectrum and this occurs for the critical 



parameter values fi c marking the threshold of linear instability. From ( 21|) such a real root 
r must satisfy 



1 + 



oo 

P.V. / dv 



r](r,/i c ) = 



oo v — r 



0. 



Corresponding to the eigenvalue A = —ikk c z is an eigenfunction 



(22) 
(23) 



iff(x,v) = e tkhcX ip(v) 



(24) 



with 



There is also an adjoint eigenfunction u) satisfying (*&, *&) = 1 given by 

Of(x,v) = ye ikkcX 4>(v) (26) 

Lj 

where 

The normalization in ( P7[ ) assumes that the root of Afc(z) is simple and is chosen so that 

< v> >= l- 

II. AMPLITUDE EQUATION ON THE UNSTABLE MANIFOLD 

Since Landau damping is weakest at long wavelengths, one expects that instability will 
occur first at k c = 2tt / L as \x is varied through \x c . This point has recently been treated 
pedagogically by Shadwick and Morrison. [34| The critical eigenvalue is then A = —ik c zo 



corresponding to Ai(z ) — for k — 1. The root z determines the phase velocity v p and 
the growth rate 7 of the linear mode 

both v p and 7 depend on /i. However it is more convenient to take 7 as the independent 
parameter and regard ^(7), 20(7), an d ^(7) as functions of the growth rate. Thus Ai(z ) = 
is understood to mean 

A 1 (*(7),„(7)) = l+ f" *!3^M> =0 . (29) 



In Appendix A, this equation is solved for 1^(7) and ^(7) to first order in 7. From (pBD , 
A = —ikcZol'j) becomes A = 7 — 2^(7) where £^(7) = kcVplj). The notation 7 — ► + for the 
weak growth rate regime always refers to the joint limit 

(7, co(j), /i( 7 )) -> (0, w(0), /i(0)) = (0, cj c , // c ). (30) 
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Subsequently, the 7 argument for v p , oj, zq, and fi will generally be suppressed when it is 
not explicitly required. 

I assume that zq is a simple root: 

Ai(*b) = ^(*b) * 0; (31) 

in addition, at z = z all derivatives of A)~(z) are assumed to have finite limits: 
lim 7 ^ + |A^(^b)| < 00 where 

A u) (z) = ^Ak (z) _ r dv vM m) 

{Z) ~ dP {Z) ~ k 2 loo dV (v - z)^ ■ {62) 
From (120|) and (p9|) Afc(z ) can be evaluated for arbitrary k 

A fe (z ) = (33) 

this identity will be needed below. Since Aj.(z) and efc(z) are identical for Im z > 0, the 
relations fl3~T|) - ([$3|) are unchanged if A k (z ) is replaced by e k (z ). 
The unstable mode corresponding to Zq is 

v) = e ik < x Mv) = e ik ' x ( -V^'VA . (34) 

V V - Zq J 

When Fq(v,h) lacks reflection symmetry, then this wave typically has a non-zero phase 
velocity and A is complex. In this case the identities A k (z) = A_k(z) and Afc(z)* = A k {z*) 
imply three additional modes: and ^ where 

M*,v) = e*" (^4). (35) 
V v-z J 

These eigenfunctions correspond to eigenvalues A*, —A*, and —A, respectively, and fill out 
the eigenvalue quartet characteristic of Hamiltonian systems. 

In the event that F (v,fi) is reflection-symmetric in v, then both real and complex 
eigenvalues may occur. If A is complex, then since k and £ commute, ^1 and 

(k-V 1 )(x,v) = e- ik < x M-v) (36) 
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are linearly independent eigenvectors for the same eigenvalue A = —ik c Zo- Thus A generically 
has multiplicity two. The same considerations hold for and ^ so the entire quartet 

has double mutiplicity. When A is real as in the symmetric two-stream instability, then 
one can show that k ■ = ^ and A is again multiplicity two. The eigenvectors ^2 and 
k . ijr 2 — ijr* correspond to —A. 



with parameter set fi = (n,u p ,Ub, A). In this example, one component, the plasma, has 
density uuq and the second component, the beam, has density (1 — n)rio; each component 
has its own drift velocity and the beam has thermal width A. The thermal width of the 
plasma has been taken as the velocity unit. If n = 0.5, A = 1 and u p = —Ub, then the family 
has reflection symmetry. 

In the four-dimensional parameter space the threshold of linear instability corresponds to 
a three-dimensional surface which is denoted by [i c . The equilibria Fo(v,fi c ) on this surface 
have eigenvalues embedded in the continuous spectrum. 

An illustrative realization of a one mode beam-plasma instability is shown in Fig. 1 for a 
system with L = 2ir, n = 0.8, u p = 0.0, A = 0.3, and Ub varied to produce the instability of 
the k = 1 mode. The spectrum of C for the stable (a), critical (b) and unstable (c) regimes 
is illustrated in Fig. 2. 



In this paper only the simplest instabilities having two unstable eigenvectors are consid- 
ered. This setting nevertheless encompasses both the case of a reflection-symmetric insta- 
bility with a real eigenvalue (two- stream) and the case of a complex conjugate eigenvalue 



A. Prototypical example 



A convenient and explicit family of equilibria, satisfying ((TTJ) , is 




(37) 



B. Critical linear modes 
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pair without reflection symmetry (beam-plasma). In either case, the components of the 
distribution function along the critical eigenvectors and are separated out by writing 

f(x, v, t) = [A(t)-9f(x, v) + cc] + S(x, v, t) (38) 

where A(t) = (\I/, /) is the mode amplitude for \l> and (\I/, S) = 0. In ( ffgp the subscript on 
has been dropped, and \I/ = exp(ik c x) ip c /L is the adjoint function for zq given in (p6|) . 
The action of the translations T a and reflection k on the distribution function implies an 
action by these operators on the mode amplitudes: from fl3lf) we have 

T a ■ A = e~ ikca A (39) 
k, ■ A = A*. (40) 

When these transformations are symmetries, these relations are useful for organizing the 
amplitude expansions below. 

The Vlasov equation (pT5|) determines the dynamics for A and S: 

A = \A+$,X(f)) (41) 

(42) 

where 

[*Mf)) = T E 7 < Wofi-i > / dv'Mv'). (43) 

L t=—oo 

In writing (§!]) I have used the adjoint relationship (^,CS) = (C^,S) = A*(^, S) = 
and in ( ^3|) an integration by parts < ij; c ,d v fk-i >= — < d v ip C) fk-i > moves the velocity 
derivative in ([L9D onto i]) c . 

C. Amplitude equation on the unstable manifold 

In ( |4"1"D and ( f42| ) the critical modes are linearly decoupled from the other degrees of 
freedom but remain coupled to S through the nonlinear terms. For 7 > 0, by restricting 
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CS + M(f)- (*,JV(/))* + cc 



to the dynamics on the unstable manifold, one can decouple the nonlinear terms as well 
and obtain from (|4"T| ) an autonomous description of the dynamics of A as a two-dimensional 
flow. This reduction to a two-dimensional submanifold is analogous to the familiar procedure 
of center manifold reduction in dissipative bifurcation theory; here the unstable manifold 
partially compensates for the absence of a low-dimensional center manifold at criticality. 

The essential properties of an unstable manifold W u are briefly described, more detail can 
be found in the extensive dynamical systems literature. |12| , |T3| The unstable modes \1/ and ty* 
span a two-dimensional unstable subspace E u and the remaining spectrum of £ determines 
the center subspace E c (for spectrum on the imaginary axis) and a two-dimensional stable 
subspace E s , spanned by the two stable modes; see Fig. 2(c). These subspaces are invariant 
under the linear flow f35| , dtf = Cf, but this invariance is lost when the nonlinear terms 
A/"(/) couple E u to E c © E s . However there are nonlinear manifolds present for the full 
nonlinear flow that are analogous to the subspaces of linear theory. Specifically, the unstable 
manifold W u is invariant under the nonlinear evolution, d t f = Cf + A/"(/), and tangent to 
E u at the equilibrium F (v, /z); hence W u is also two-dimensional. Solutions on this manifold 
f u (x,v,t) asymptotically approach F as t — > — oo. 

The invariance of the unstable manifold means that restricting the Vlasov equation to 
W u yields an autonomous two-dimensional dynamical system describing the evolution of 
initial conditions on the manifold, f u (x,v, 0) G W u . Near the equilibrium this restriction is 
tractable since the tangency between W u and the unstable subspace allows the manifold to 
be described by a function, 

H :E U -> E C ®E S (44) 
(A, A*) -> H(x,v,A,A*) (45) 

which measures the "distance" from E u to W u ; see Fig. 3. Thus near Fq, for trajectories 
f u (x,v,t) on W u , the evolution of the non-critical modes S(x,v,t) in (|3"5|) is controlled by 
the critical modes: 
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S(x, v, t) = H(x, v, A(t), A*(t)), 



(46) 



and these trajectories can be described entirely in terms of H and the evolution of A(t): 



f u (x, v, t) = [A(t)*(x, v) + cc] + H(x, v, A(t), A*(t)). 



(47) 

If H is known, then using (|47"D , the general equations fl4~i~D and (|42"D can be restricted to 
the unstable manifold: 



dS_ 
dt 



A = XA + (^f,J\f(f w )) 

= £H + Af(f u ) - f Af(f u )) * + cc 



(48) 
(49) 



Now (|48D defines an autonomous two-dimensional flow describing the self-consistent nonlin- 
ear evolution of the unstable mode; this is the amplitude equation I wish to study. 
Translation symmetry forces the right hand side of fl48|) to have the form 

XA + ^,Af(r))=Ap(a^) (50) 

where a = \A\ 2 and the function fx) is not constrained by the translation symmetry. |55fl 
Typically p(cr, /x) is complex- valued, however when Fq is reflection-symmetric then p(a, /i) is 
forced to be real. 

It is convenient to view (]50f) as expressing p in terms of H, and make this connection 
more explicit by evaluating (ty,J\f(f u )). This calculation requires the Fourier series for H 



H(x,v,A,A*)= J2 e ikk * x H k (v,A,A% 

k=—oo 



and the Fourier components of /" 



= [A1> c (v) 5 k>1 + A*ip c (v)* + H k (v, A, A*). 



(51) 



(52) 



The evaluation of (ty ,J\f(f u )) from (H3I) is simplified by first noting that the components of 
H are forced by translation symmetry to have the form 



H (v,A,A*) = ah (v,a) 

Hi(v, A, A*) = Aahi(v,a) 

H k (v, A, A*) = A k h k (v,a) for k > 2, 



(53) 
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where H_k = H k and the functions h k are not constrained by the translations. [[37| However, 
if reflection symmetry also holds, then 



Combining (|52] 



h k (-v,a) = h k (v,a)*. 



with (^,7V(/ U )) in (0) yields 



— — < < d v ip c , (h - h 2 ) > +y < d v $ c , i>* c > 



(54) 



(55) 



-a 



00 ^-2 



Z=3 



where, on the right hand side, the velocity integral of h k is denoted by 



T k (a) = 



dv' hk(v', a). 



(56) 



Now comparing (|50|) and ([55D provides the desired expression for p 
p(a, fj) = X - y < < d v ip c , (ho - h 2 ) > +y < d v ip c , ip* > 



(57) 



r\ < <9„^ c , /i > -r; < <9„^ c , /i 2 > +y < <%^ c , ^ > — ^ < d v ip c , h 3 > 



00 1-2 



+ E ^i r « < <^ c , K-x > -vr; < d v ^ c , h l+1 >] 



1=3 



In order to exploit this expression for p, it is necessary to determine H or equivalently to 
determine the functions h^. 

D. Representation of the unstable manifold 



An equation for H follows by requiring consistency between (f46|) and (fK|); setting the 
time derivative of fllB]) equal to the right hand side of (|||) gives 

^_A+—A*=£H + N(f u ) - [(*, N(f u )) * + cc] (58) 

which is to be solved for if subject to H (x, v, 0, 0) = and 
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dH , dH 



(59) 



These latter conditions are implied by the tangency between W u and E u at the equilibrium 
(A, A*) = (0, 0), and are automatically satisfied in this case by virtue of (|53|). 

With the previous expression for (ty,J\f(f u )) in (^) and the notation in (^) for the 
Fourier components of H, the components of (EW) take the form 



dA OA* 
ia d 

k c dv 



CO l-l 

r c + <r{h\ - 4> C T\) + a 2 h*^ + E — fyT, 

1=2 ' 



CC 



(60) 



CM 



<9A* 



L4<t d 
k r ± dv 



h -h 2 + -ip* r 2 + a 



hoF, - h 2 T\ + -h\T 2 - -h 3 T* 2 



(61) 



oo 1-2 



1=3 



CM 



CM* 



ZA 2 9 
fen <9t> 



hi + ^cTi -h 3 + ^h T 2 + 7^c r 3 



(62) 



and 



a' 

v 3 



6 1=4 1 



dH k 
dA 



A 



dA* 



A* — L k H k 



%A k d . , 
-T- ^- { h k -i + 



k-i 



k ~ 2 h 

r*-i + E 

Z=2 



"fc-l-L 1 — n k+l + 1 7 1 fe-1 + I -1 + ~, — —T' 

K — 1 K K + 1 



fe+1 



-h k+1 ri + 



k+ 1 



fc+i 



+ E 



oo (jl—k 



l=k+2 



(63) 



oo ^ 



E -rhk+iVi 



1=2 
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for k = 0, 1,2, and k > 2, respectively. In the k = 1 component (|61|), P is the projection 
operator onto the ip c (v) component of a function g{y ), 



(64) 



and the orthogonal projection is denoted by V± = I — V. 

The expressions for AdAH k + A*dA*H k can also be evaluated in terms of the functions 
hk{v, a) and p(a, /i): 

dh 



-3A A + m A ' = °<* + p' 



h + a 



da 



(65) 



{(2p + P *) + (p + p> 



do" 



and 



(66) 



^A+^A*-L k H k = A k 
OA OA* k k 



{kp - L k }h k + (p + p*)a 



dh k 
da 



(67) 



for k = 0,1, and k > 2, respectively. By combining (|6DD - fl6"3|) and (|65| ) - (|67]), the component 
equations of (|58| ) reduce to a simpler form: 



da 



(68) 



i d 
k c dv 



oo ^-i 



r c + <r{h\ - ^ C V\) + a 2 htTx + ]T —h^T 



/=2 



CC , 



{(2p + p*)-L 1 }/ il + (p+p>-^ 



^4 k-^ + ^:r 2+ . 



(69) 



/nA - h 2 v\ + ^r 2 - h 3 r* 2 



00 a 1 ' 2 r 

+ E-r-K-i^-^ + irr 



{2p-L 2 }/i 2 + (p + p> 



9a 



(70) 
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i d 



h x + ip c r 1 -h 3 + i/i r 2 + ^* r 3 



a' 

r3 



a" 00 rr l ~ 2 r 

d (=4 4 



and 



{/cp-L fc }/i fc + (p + p*)o- 



dhk 
da 



k - 2 h 



1=2 



-a 



hk-iTi — hk+i 



Ifc-i + — Ifc + - — — -Tk+i 



k-1 



k 



-hk+iT 1 + 



k+1 



fc+i 



00 Q-l—k 



E 



k+ 1 



l=k+2 



00 ^ 



1=2 



(71) 



for A; = 0, 1, 2, and k > 2, respectively. 

Together with (|57| ) these component equations determine the functions p(a, fi) and 
/i)}£L . From a practical point of view, all that has been achieved to this point 
is a reduction of the problem to the analysis of functions of a single real variable, i.e. a. 



However, in the study of the amplitude equation p8|), this reduction does provide a useful 
simplification which is exploited in the discussion below. 



III. EXPANSIONS, RECURSION RELATIONS, AND PINCHING 

SINGULARITIES 

The amplitude equation on the unstable manifold, 



A = Ap(a, /i), 



(72) 



is analyzed by expressing p as a power series in the mode amplitude, 
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j=0 

whose coefficients Pj{p) are calculated from fl5T|) using an analogous series for hk 



(73) 



h k (v,a) = h k,j(v) a 3 . 
For notation, denote the integral over hkj(v) by 

/oo 
dv h ktj (v) 
-oo 

so that 



(74) 



(75) 



from (|56|), then expanding the expression for p in (W7f ) gives the coefficients Pj(fi) as 



(76) 



A 



Pj = -j-[Aj + Bj] for i>! 



(77) 
(78) 



where 



Aj = - < d v ip c , (hoj-x - hzj-i) > -- < d v il) c , i/j* > V 2 j-i 



i~ 2 

~ {< d vi>c, h ,j-i-2 > Fiji- < d v ^ c , haj-i-2 > 

1=0 

i-* f i 
2 



J] - < ^~ c , /i*^_ 2 > r 2i/ + £ 



=0 



m=0 



— 2-fl,m . rv i r * 

7T^ < ^c, hj_ hl . 



J-1 + 



> 



"P* 

j—l 7 m 



(79) 



(80) 



Here and below, a summation is understood to be omitted if the lower limit exceeds the 



upper limit. The organization of terms between (I79J) and (50) will turn out to distinguish 



different singular behaviors in the 7 — > + limit: A,- ~ ^/ _ ( 4 J _1 ) and ~ <-y-( 4 J _2 ); thus the 
terms are sub-dominant in the weak growth rate regime. 

The leading term in ([7^) is the eigenvalue po = A, and higher order terms are determined 
by calculating hk,j(v) from (|68[) - (|71); the results are summarized below. For k 7^ the 
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coefficients hkj(v) are expressed in terms of the resolvent operator R k (w) = {w — L k ) 
an arbitrary complex number w, the resolvent acts on a function g(v) by |3£ 



(R k (w) g)(v) 



ikk c (v — iw/kk c ) 



9{v) 



rj(v,n) 



k 2 A k (iw / kk c ) 



dv' 



v' — iw/kk c 



l . For 



(81) 



A. Series coefficients for h k (v,a) 



For k = the coefficients h Q j are found from flog) ; inserting the expansions for h k (v,cr) 
and rfc(cr) into ( j6S|) and setting the coefficient of cr z to zero yields 



hoAv) 



(1 + /)(A + A*)' 



The functions Iqj(v) are given by 



k c dv 



and for / > 1 



z-i 



Jo,j(«) = " E (! +i)(Pi-i + M u ) 



i d 
k c dv 



1-2 



3=0 



l-l 3 h* 

^ ^ I — i + 1 



cc 



(82) 



(83) 



(84) 



For k > 1, following the same procedure in (|HD| ) - ( |7T|) determines the corresponding 
expressions for /i^/; these coefficients have the form 



h k i(v) = R k {w kt i) I kl 



(85) 



where R k (w) is given in (plf ) and 



w fci/ = (Z + 4,i) (A + A*) + k\ = 2(1 + 4,i) 7 + feA. 



(86) 



For k = 1, the functions Ii,i{v) are given by 
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h,i(v) = -E [(2 + j)m-i + (1 + j)pl 

3=0 

+f V±4- \h ,i - h 2 ,i + \rj2,i 



k r dv 



l-l 

+ E 

3=0 



ho,j^i,i-j-i - h 2 ,jT* 1 i_ j _ 1 + -hl tj T 2 ,i-j-i 



-l^ 1 ^ 2,1-j-l + E 



2, //l 



m=0 



7— j+l,m ■ 

J-J + 2" 



l—j+2,j-m 



1-2 j 

EE 

j=0 m=0 



hl-j+2,m p» 

Z-J + l ^' +1 J- m 



For k = 2, from (|70|), the functions I 2 ,i{v) are given by 



for / = 0, and for / > 1 
l-l 



hiiv) 



E V- • ./)/'/ , • .//'/' , 

3=0 



L 2,j 



l-l 



id] 11 

+ F dv V 11 ' 1 ' 1 + ^ cTl ' 1 - 1 ~ hs ' l ~ 1 + 3^* r 3,«-i + - E h o,j T V-j-i 



1-2 



1 1 

+ E h\jY\).-j-2 - h^* u _ 3 _ 2 - - h ^h-j-2 + q^i,i r 3,/- 

J=0 



+ E 



2 



l-j,m 



' l-j + * 



l-j+2,j-m 



1-3 h 1-4 j h 

ll.1l'.,. x x Ill-iA 



i=0 u j=0m=0 ^ 

For k > 2, the functions Ik,i(v) are given by 

4,i0) = - E i( k + 3)Pl-3 + J Pi-j] h k,j(v) 



3=0 



+ 



z d 
k r dv 



■0, 



fc-2 I 



hk-i,i + t -Fk-i,i + E E 

fc ~ 1 l'=2 i=o 



hk-v 



—hk+i,i-i + 7 — — rfe + i;_i. 
K + 1 

+ E 

J=0 



k,i-j 
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1-2 



+ E 

3=0 
3 

+ £ 



( K 



h* 



r 



-Jjn 



r 



l+k—j,j—m 



k+l,l-j-2 



hk+l—j,m -p* 

l-j l -^-m 



in this last expression, if a subscript is negative the term is understood to be omitted, e.g. 
for I = 0, hk+i,i-i is omitted. 

The arguments w k ,i to the resolvent fl8"5|) determine poles of hk,i(v) located at v — z k j 
where z k ,i = iw kj i/kk c - For > 1, these poles always fall in the upper half plane above the 
phase velocity v p : 



z k,l — z + 



ijd, 



k.l 



V v + 



17(1 + d ki ) 



(91) 



k c k c 

where d k j = 2(1 + 8 k> i)/k. 

The relations ( ff8| ) and ([82]) - fl9~0D can be solved systematically to calculate pj and ft,*^ to 
any order. The leading coefficient po is determined by linear theory (\T7\) and from the linear 
eigenfunction ip c one can also calculate /i ,o an d h 2> o, c.f. ( p3|) and (|88|), respectively. These 
two coefficients then suffice to calculate p\ from (f78|). From {pi, /i ,o, ^2,0}, the coefficients 
/i^o and can be determined and then h 01 and h 2t i- This provides the input to calculate 
P2, and from this point on the structure of the calculation to all orders falls into a recognizable 
pattern that is summarized in Table I. 

The expressions (|82| ) and (|85|) for h kj i allow a useful evaluation of T k i in terms of I kj i. 
For k = the coefficient vanishes identically To,; = 0, and for k > 0, from (|75|) and ([So]), 
one has 



dvR k (w k ,i) I, 



k.l 



dv 



T)(V,H) 



k 2 A k (z k ,i) J-oo v - z k ,i 
k 2 A h (z k>l ) - (Ai(z k>l ) - 1) 



k 2 A k (z k>l 



-oo f' - z fc ,« 
, h,i(v') 



dv 



V - z k i 



00 dv' W) 



V - z k . 



(92) 
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where in the last step the identity k 2 Ak(zk t i) = k 2 — 1 + A.i(zk,i) has been used. The final 
expression in ([92] ) will prove helpful in analyzing the form of I\z as 7 — > + . 



B. Analysis of the cubic coefficient 

It is instructive at this point to evaluate and examine the cubic coefficient in detail. 
From fl78[) we have 



< <9„^ c , (/i 0j0 - h 2 ,o) > +7; < dvipc,4>c > r 2 , 



(93) 



so that ho,o and h 2 fl must be found from flS2"|) and (j55|). A straightforward calculation yields 



_(v-zq) (y-z%) 



(94) 
(95) 



(v — Zq) 6(v-z ) 

where A^(zo) is the second derivative defined in (|32|). Integrating these functions over 
velocity yields r ,o = and 

-AS 2) (^o) 



r2.o 



3 k 2 . 



(96) 



The remaining integrals in ( p3|) are expressible in terms of the derivatives of Ax(zq): 



< d v ip c , if>* > = — 
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< d v ip c , hop > 

< d v ip c , h 2 ,o > 
Thus we find from fl9"3]) 



4 7 3 



2(A?\z )) 2 -3A?>(z ) 
48A£Ai(*b) 



z 7 AS 2) (^o) 2 7 2 aP ) (^ ) 

( 4 )/ 



(97) 
(98) 

(99) 



1 



Pi 



4^3 



27Ai 2) (^o) 7 2 / 2A^(^ )-A;(^o)A^(z ) 
fccAHzo) 3^ 1 Ai(2b) 



(3), 



(2), 



17' 



Ylkl M 



^(3AS 4) (.o)-2(aS 2) (, )) ; 



(100) 
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which simplifies to 



Pi 



- 5 + 0(7) 



(101) 



as 7 — > + . 

The cubic coefficient is obviously singular in the weak growth rate regime and the strength 
of the singularity is set by the contribution from h Q0 (v) which is the nonlinear correc- 



tion to the equilibrium Fq at this order. [39j The following summary of the calculation of 
< d v ip c , ^0,0 > in ( HP pinpoints the origin of the singularity. From ( |94|) and the definition 

of ip c , one has 



1 



00 



1 







00 



-7){v,fi) 



(V - Z Q ) (V - zg) 



G?t> 



r/(w,/i) 



(102) 



fc c 2 Ai(z ) ^-°o " (v - zoY (v - 4) ' 
Since zq = v p + ry/A; c in ( |102| ), there is clearly a pinching singularity at v p when 7 — > + . 
By expanding the integrand in partial fractions, 



1 



(fc c /2* 7 ) (k c /2i 7 ) 2 , (fc c /2i 7 ) 



+ 



(v - z ) 4 (v - *3) (v - ^o) 4 (« - ^o) 3 (v - z ) 2 

4 



(103) 



1 



U - Zq V- Zl 



( |102|) becomes 

< <9„^ c , /i ,o > 



*c 2 A'i(*b) 



+ 



kc_\ 



(Ax(«b) - 1) - (Ai(^o) - 1) 



(104) 



A'i(^o) 



2i 7 



' k c \ A?\z ) 
2t7 J 



+ 



6 



Since Ai(#o) — and Ai(zq) = 0, the 7 4 terms vanish identically and do not contribute to 
the 7 — > + limit; thus ( |104| ) reduces to (|98|) and a 7~ 3 singularity. 
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The interpretation of the 7 3 singularity is suggested by examining the form of the mode 
equation (|72|), truncated at cubic order, in the 7 — > + limit: 



A = A 



(105) 



where A = 7 — iuj; in terms of amplitude and phase variables A = pe %e this reads 

1 



P = P 



6 = u-Re 



(106) 
(107) 



Ai(*b) j 4fc c7 2 

The manifest singularities at 7 = can be removed by introducing a rescaled amplitude 
variable: 

p{t) = 7 2 rirft) (108) 
which evolves on the time scale r = r yt. In terms of r(r), the mode equation becomes 

(109) 
(110) 











= r 


[i-i 


dr 




4 




= UJ 


7 2 




4fc c 



aS 2) (^o) 



Ai(^o) 



r 2 + 0( 7 V). 



For the terms shown explicitly in ( |109|) - ( |110| ), the 7 — >• + limit is now well behaved; this 
suggests that the effect of the singularities in the coefficients is to produce the scaling p ~ 7 2 
in the dynamics of the mode amplitude. Our ansatz in (|108|) bears some resemblence to the 
simpler case of Hopf bifurcation discussed in the Introduction but the scaling is quantitatively 
different indicating a much stronger nonlinear effect and a smaller nonlinearly evolving mode. 

It is also important to note that to this order the result ( |108| ) has a remarkable degree of 
universality. The 7~ 3 singularity in ( |101| ) sets the scaling in ( |108j ) and the coefficient, —1/4, 
of this singularity is completely independent of the underlying equilibrium F (v,p c ). Thus, 
for example, the singularity is the same for a beam-plasma instablility (complex A) as for 
a two-stream instability (real A). Indeed to this order, the rescaled amplitude equation for 
dr/dr is also completely independent of F as 7 —>■ + . 
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These remarks tacitly assume that the neglected higher order terms in ( |105|) do not alter 
conclusions reached for the truncated equations. The remainder of the paper is devoted to 
a systematic analysis of the higher order terms in the series for p and hk,i- It will turn out 
that the coefficients pj for j > 2 are also singular and that the naive estimate 0(7 8 ) in 
( |109| ) and ( |11(J| ) is not correct. In fact the terms that are higher order in r are not higher 



order in 7, rather they appear to be as important for the dynamics of A(t) as the cubic term 
A|A| 2 . Nevertheless the singularities to all orders are absorbed by the scaling in ( |108| ). The 
dependence of the higher order terms on F (v,fi c ) as 7 — > + will be discussed later. 

IV. SINGULARITY STRUCTURE OF THE EXPANSION 

The detailed calculation of pj rapidly becomes prohibitively laborious, but the recursion 
relations determining the higher order coefficients in terms of lower order quantities (cf. 
Table I) can be analyzed to obtain useful information. Most basic is the question: for j > 1, 
how singular is pj as 7 — > + ? This issue requires an accurate estimate of the build up of 
pinching singularities in the integrals found in pj. For this purpose I introduce an "index" 
which allows the divergence of a given integral to be assessed by a simple counting procedure. 

A. Definition of the index 

For n > define 

D n (a,v) = - - — r (111) 

(v - ai){v - a 2 ) ■ • ■ {v - a n ) 

where a = (ai, . . . , a n ) and define D (a, v) = 1. Evaluating pj for j > 1 involves integrands 
of the form 

Q(v, /i) = DM v)* D n (a, v) ^ ^ (112) 



with m + n > 1. The poles in ( |112| ) are given by 
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acj = zq + i^Uj/kc j = 1, ... ,n (113) 

/3* = z* - z 7 & IK j = l,...,m; (114) 

hence they lie along the vertical line Re v — v p at locations determined by the numbers 
Uj > and Q > which are assumed to be independent of Fq for all j; in particular Uj and 
Q are independent of 7. The function Q(v,[i) depends on /1 through zq, 7, and n(v,fi). 
The index of Q(v,fi) in (|112|) is defined to be 

Ind [(?] =m + n + /-2. (115) 

Since we assume m + n > 1, Ind [Q] > — 1. If mu 7^ 0, then as 7 — > + , the integral of Q 
diverges due to a pinching singularity at the phase velocity v p . When Ind [Q] > 0, the index 
of Q gives the maximum possible strength of this divergence. 



Lemma IV. 1 For Q(v,n) in ( 112 ) with m + n > 2 and J = Ind [Q], the integral of Q 
satisfies 



lim 

7^0+ 



7 J / dvQ(v,fj.) 



oc 



< 00. (116) 



// J replaced by J — 1, then the limit ( |116|) diverges in general unless mn = in which 
case the limit is zero for any J > 0. 

Proof. 

Since 

|" (foA»0M*A»(a,tO = (117) 

/■oo $- X T)(v ll) { m n 

i-00 dV dv 1 - 1 VpD^&favY D n (a,v) + J2D m ((3,v)* D n+1 (a, aj ,v) 

we can reduce to the case with 1 = 0. If mn = then there is no pinching singu- 
larity and the integral in ( |116| ) has a finite limit, cf. (0); hence the limit is zero if 
J > 0. Assume that mn > 0, then for integrals Jf^, dv D m (j3, v)* D n (a, v) n(v, fi) 
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with m + n = 2 the integration can be easily done and ( |116| ) explicitly verified. 
For integrals with m + n > 2, by expanding the integrand in partial fractions, 
they can be re-expressed in terms of integrals with m + n — 1. A simple induction 
argument then establishes ( |116| ). In Appendix B the limit in ( |116| ) for I = is 



shown to be non-zero in general so, barring an accidental cancellation, if J is 
replaced by J — 1 then the modified limit will diverge as 7 _1 . □ 

The definition in ( |1 1 5| ) can be generalized in some obvious ways without sacrificing its 
usefulness. Suppose q(/j,) is a function of \x with the asymptotic behavior 

lim g(/i)~ 7 " 5 , (118) 

7^0+ 

then we define the index of q(fi) Q(v, fi) to be 

Ind [qG] =Ind [Q] +5. (119) 

Estimates of the form ( |116[ ) still hold for qQ with J = max {6, Ind [<?{?]}; the anomalous 
case J = 5 arises when Ind [Q] = —1 since the integral of Q is then nonsingular and the 
asymptotic behavior of q J dvQ is determined by g(/x). Finally if Gi(v,p,) and ^2(^5 have 
indices satisfying Ind [Qi] > Ind [Q 2 ] then we define the index of the sum to be the larger 
index: 

Ind [Q x + g 2 ] =Ind [&]. (120) 
B. Examples 

As a simple example of the index, note that the functions ip c , h^^, and /i2,o have indices 
given by 

Ind fok] = -1 (121) 
Ind [ho tQ ] = 1 (122) 
Ind [h 2>0 ] = 1; (123) 
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similarly the integrands in fl97|) - (0) give 

Ind [d v ^*^ c \ = 1 (124) 
Ind [d v $*ho t0 ] = 3 (125) 
Ind [d v <$*h2 )Q ] = 3. (126) 

From Lemma IV. 1, this information immediately tells us that the singularity of pi cannot 
be worse than 7~ 3 . 

A further useful observation is the effect of various operators on the index. 

Lemma IV. 2 For Q(v, fi) as in (|112|) , the operators d v G, V±G and Rk{wk,i)G either leave 
the index unchanged or else increase it by one: 

Ind [d v G] = Ind [Q] + 1 (127) 
Ind [V±G] = Ind [Q] (128) 
Ind [R k {w k>l )g\ = Ind \g\ + 1. (129) 



Proof. 



The first identity is obvious, and the second follows immediately from the defi- 
nition of V± below (|B3). For the third identity, the resolvent (|5TD gives 



Rk{wk,i) Q 



1 



ikk c (v - z k j) 



k 2 A k (z k: i) J-oo v' - z k: i 



dv 



(130) 



For the first term the index is increased by one, Ind [G/(v — Zk,i)] — bad [G] + 1. 
In the second term, Ind [rj(v, fi)/(v — z k j)] = — 1, so the function multiplying the 
integral 



?0) 



1 



k 2 k k (z k ^ 



°° dv' Q -^ 

-oo V — Z k I 



(131) 



will determine the index and the singular behavior of the corresponding velocity 



integral (recall the discussion below (|119|) ). Since 
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-CO f ' — Zk t l ry 



1 



(i+Ind [Q]) 



and 



k 2 A k (z k>l ) = k 2 - 1 + Ai(^,,) = /c 2 - 1 + A; (,8b) 



(132) 



(133) 



g(/i) has the asymptotic form 



lim 

7^0+ 



7 J q(») 



< oo 



(134) 



with J = Ind [Q] + 1 + 5fc 5 i. Hence the index of the second term in (|130 ) is 
Ind [Q] + 1 and Ind [Q] for k = 1 and k ^ 1, respectively. In either case, the 
overall index is given by ( |129|) . 

□ 



C. Estimate of singularities 

The main result on the singularity structure of the amplitude expansions for p(a, p) and 
h k) i(v, cr) can now be proved. 

Theorem IV. 1 For k > and I > 0, the indices of I k i and h k i are given by 



Ind = Ind [h ktl ] - 1 

Ind [h Kl ] = 2k + 41 - 3 + 4(5 fc ,o + M- 



(135) 
(136) 



The asymptotic behavior of the coefficients F k j and pj satisfies 



lim 

7^0+ 



7 J T k)l < oo where J = Ind [h ki i] + 4,1 



lim 

7^0+ 



7 4 - 7 1 Pj(f-i) < oo for j > 1. 



(137) 
(138) 



Proof. 
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1. Since ho,i = Io } i/2j(l + I) and h k ,i = Rk(wk t i)h,i, it follows immediately 
from ( |119| ) and ( |129| ) that the indices of and hj-j differ by one as stated 



in ( |135| ) - provided the index of I k ,i is well defined. The induction argument 
outlined below shows that each function I^i is obtained via the recursion 
relations as a sum of terms each having a well defined index. Thus 1^ will 
always have a well defined index. 



2. The estimate in (|137|) follows from the formula for Tk,i in (|9~2"D: 



k l - 1 + Ai(z fc) j) J-oo v-z k ,i 



° dv M±„L } (140 ) 



as 7 — ► + , the integral cannot diverge more strongly than 

f Ind [/i fc ,i] 

-oo v - Z k>i \7 

and the function multiplying the integral is nonsingular unless k = 1 in 
which case it is (9(7 _1 ). 

3. The index formula for h k ,i ( |136|) and the estimate for in (|138|) are 



proved by induction using the recursion relations. The induction argument 
is organized by the pattern shown in Table I. At the top of Table I, the 
results in ( |136| ) and ( |138| ) have been explicitly verified for {/io,o> ^2,0} an d 
Pi in Section III.B. 

4. Moving downward in Table I, the next quantities to consider are h\fl and 
hzfl from (|87|) and (^3), respectively: 

hi, = ^iKq) {/i ,o - ^2,0 + ^:r 2 ,o} (141) 

^3,0 = ^Ko) |j 1^2,0 + yr 2l0 | • (142) 

Reading right to left in these expressions, the index for each of the quantities 
in braces {• • •} is 1, and this index is increased to 2 by d v and further 
increased to 3 by the resolvent; hence Ind [/i^o] = 3 and Ind [h^o] = 3 in 
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agreement with ( |136| ). Next, according to Table I, are ho,i and h 2j i, from 
(B3) and PP they are: 



+ 



-(pi+p*)/io,o(u 
z 9 



(143) 



ft* 



— cc 



ft-2,1 = R2(w 2 ,l) 



-2p 1 h 



2,0 



(144) 



z 9 f 1 I l 



For /io,i the expression in braces {• ■ •} has index 3 and the larger expression 
in brackets [• ■ ■] has index 4, so Ind [/io,i] = 5 in agreement with ( |136| ). 
Similarly the index of h 2t i is seen to have the correct value Ind [h 2t i] = 5. 

5. This provides enough information to verify ( |138| ) for p 2 . From (^) this 
coefficient is p 2 = i[A 2 + B 2 ]/k c where 



A, 



Bo 



< d v tp c , (/i ,i - h 2A ) > +^ < d v 4} c ,ip* > r 2 ,i 



+ < d v ip c , h , > Tifi- < d v ip c , h 2 ,o > 



< d v ip c , h* lfi > r 2i0 + < d v ip c , h* 2fi > 



(145) 



(146) 



1 2,0 



< d v ip c , h 3 ,o > 



Since multiplying by d v ip c increases the index of an integrand by two, cf. 
( |121|) - ( |126| ), it is straightforward to verify that p 2 satisfies the estimate 
in ( |138| ). For example, consider the first term in (|145|) . Since Ind [/io,i] — 
Ind [/i2,i] = 5, the integrand has index Ind [d v ip c (/io,i — ^2,1)] = 7 so the 
worst possible singularity for the integral is 



< d v ip c , (h ,i - h 2 ,i) >' 



r 



(147) 



which is consistent with ( |138| ). A similar verification for the remaining 
terms in p 2 is straightforward. 
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6. According to Table I, for fixed N > 1, given the functions 

{h k>l (v) | k = 0, . . . , N+ 1 and I = 0, . . . , l max {N, k)} (148) 

where l ma x(N, k) = N — k + 1 — 2(5^0 + <5fc,i) and the lower order coefficients 

{ Pj \j = l,...,N}, (149) 

then one can calculate the additional functions 

{h k>l {v) | k = 0, . . . , N + 2 for I = l max {N + 1, k)} (150) 

and the next coefficient Pn+i- For N = 2, it has been shown that the func- 
tions ( |148| ) and the coefficients ( |149| ) satisfy ( |136| ) and ( |138| ), respectively. 



7. The assumption for the induction step is now clear. For arbitrary N > 2, 
assume that the functions (|148|) satisfy ( |136|) and also that the coefficients 



( |149|) satisfy the estimate in ( |138j ). Then from this hypothesis, the recursion 
relations determining the functions in ( |150|) and for Pn+i imply that these 
new quantities will also satisfy (|136|) and ( |138j ), respectively. Verifying this 
conclusion involves only a routine evaluation of the indices arising from the 
recursion relations as was done above for iV = 2. This completes the proof 
of O and (|T38|) . 

□ 

D. Implications of Theorem IV. 1 

The index formula for h k ^{v) and the asymptotic estimate on T^/ imply that in the 
the general expression for pj (|78|) the dominant terms are in Aj. More precisely, one finds 
Aj ~ ^~( 4 J -1 ) and Bj ~ ^~( 4 J~ 2 ) as 7 — >• + . The estimates in (|137| ) and ( |138|) motivate the 
following asymptotic definitions for I\z and pf 
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lim 

7~>0+ 



7 



(2fc+4«-3+5<5fe,i) r 

1 k.l 



bj{nc) = lim [ 7 4j 1 p j (fi) 



lim 

7^0+ 



for k > 1 and / > 



7 4j '- 1 Aj(jJ.j\ for j > 1. 



In this notation, the asymptotic form of for k > 1 is 



I 

,7, 



+ Oil)} 



where 



(151) 
(152) 



(153) 



1=0 



21 



(154) 



1. Mode amplitude equation 



The full mode amplitude equation, 



d A 00 

— = XA + AY,PiM° i (155) 

when written in polar variables A = p e~ %e gives 

do 00 

f ^pj + P Y[Re Pj (n)]p^ (156) 
-^^-^[Im^p^ (157) 

The 7 — > + limit can be taken using the asymptotic behavior for pj (|152|) given by Theorem 
IV. 1 and the rescaling p{t) = •y 2 r(jt) suggested by the low order analysis (|108|) : 



1 + E [Re6,(/i c ) + 0( 7 )]r 2 ^ (158) 
3=1 



dr 
dr 

df) 00 

a - = u - 7 £ [Im &^ c ) + 0( 7 )] r 2 ^ (159) 

where = fc c f p is the linear mode frequency. These equations are non-singular at every 
order as 7 — > + , and the final asymptotic equations are 

dv 

— = rR(r,p c ) (160) 
d6 , , 

(161) 
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where 



R(r, t M e ) = {l + 'ir i \Reb j (p e )] 



(162) 



3=1 



If the equilibrium is reflection-symmetric, then u = and [Im bj(fi c ) + 0(7)] = so the 
phase equation ( |159[ ) becomes d6/dt = 0. 

Note that, in contrast to the rescaled amplitude equation for Hopf bifurcation (^), in this 
case the higher order terms in r are not higher order in 7; thus there is no apparent justifi- 
cation for a truncation of the expansion after the leading nonlinear terms. The dependence 
of R(r, /i c ) on /i c is analyzed in Section V. 



2. Electric field 



The electric field, 



E(x,t) = Y,(e tkkcX E k (t)+ cc), 

k=l 



(163) 



is obtained from the potential E = — d x &, and from Poisson's equation (|13|) the Fourier 
components are 



ikk c E k (t) 



dvft(v,t). 



(164) 



With fl52|) for f%(v,t) this becomes 



E k (t) = I 



-A(l + aF 1 ) k = 1 



(165) 



hence the field is given by 



E(*, t )=[ T 



00 _,2fc-p „fc 



7 2 r ( r ) (1 + 7 4 r 4 ri ) e i(fc c x-e(t)) + 1 k c ik(k c x-6(t)) 



k=2 



+ CC 



(166) 



In the 7 — > + limit, this gives 
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E(x,t) 



r(r) l + rjr(r) 4 + 0( 7 ) e 



D i(k c x-6(t)) 



(167) 



7E 



2E. (T c k r(r) k 



0( 7 ) 



D ik(k c x-e(t)) 



k=2 



k 



+ cc 



where T£ is defined in ( |154| ). This result describes an electric field that obeys trapping 
scaling E ~ 7 2 and also indicates that all higher spatial harmonics are uniformly 0(7) 
relative to the wavenumber of the unstable mode. 



V. ASYMPTOTIC RECURSION RELATIONS 

From the definition of bj ( |152| ) one expects these limits to depend on the parameters fi c 
through the underlying equilibrium F (v,fi c ). A striking feature of the explicitly calculated 
cubic coefficient ( p. 01] ) was the result 

h = -i; (168) 

thus bi is a pure number and independent of fi c . It is natural to ask what happens for the 
higher order singularities. The dependence of bj on \l c can be systematically investigated by 
taking the 7 — > + limit of the recursion relations directly rather than calculating individual 
coefficients pj and studying their asymptotic behavior. 

Let hk,i and Ik,i denote the terms of maximum index in hkj and Ik,i respectively. From 
Theorem IV. 1, this definition implies 

hk,i = hk,i + [terms of index 2k + 41 — 4 + 4(5^0 + 5fc,i) or less] (169) 
h,i = h,i + [terms of index 2k + 4/ — 5 + 4(5^0 + 5fc,i) or less] (170) 

Since only terms of maximum index can contribute to the limit in ( |152|) , the expression for 
Pj in (|78| ) yields an explicit formula for bj 



ik c bj = lim 7 4j 1 

7^0+ 



< d v ip c , (hoj-i - > (171) 

+ Yl (< dvipc, hoj-1-2 > r^- < d v i/; c , h 2 ,j-i-2 > r*^ 
z=o 
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Thus to study bj, only recursion relations for hk,i are required. 



A. Truncated recursion relations 



From (55) - (p0|), these relations are obtained by simply discarding terms of submaximal 
index. This truncation yields the following expressions: 

hi 



h ,i 



(1 + 0(A + A*) 



(172) 



for k = 0, with 



hi( v ) = < 



I = 



(173) 



- rfi (1 + j)(pi-j + PUj) hoj(v) I > 

+1 1 { (ki-x ~ Wi-i + E;= 2 o ^IW-a) - cc} 
For k > 1, we obtain hk,i{v) by calculating Rk{wk,i) h,i and discarding any terms of index 
less than Ind [hk,i\- The required forms for If.j are 



/-i 



id' 



^2,jT* 



j x M-i-i 



(174) 



-E[(2+j>-,- + (i+iM-i] k 



3=0 



hi(v) = < 



1 = 



h 2 j I > 



(175) 



+il{W-i-Vi + ^riM 

+ Ej=0 ^1,5^1,1-^-2 - ^3,3'r^_ j _ 2 



} 



and for k > 3 



4,z(v) 



z d 
k r dv 



l-l ^ 1-2 

h>k-i,i — hk+ij-i + Z ^k-ij^iji-j-i — Z hk+i,jFlj_j_ 

j=0 3=0 
l-l 

3=0 



(176) 
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B. Integrated recursion relations 



The expression for bj ( |171| ) depends only on certain integrals of h kt i which can be obtained 
directly from an appropriate integrated form of the truncated relations. For non-negative 
integers (m, n), consider limits of the form: 



kj 7-0+ ^ A; (z ) 

m+n— 1 / J 



^ v m+n-2 / 

S(m,n;P,a) = [ — ) lim — — — / dv D m (P, v)* D n (a, v ) rj(v, jj) I (177) 



B k ,i(m,n;P,a) = — lim / du D m (/3, «)*£)„(«, u) /i fc l (u) (178) 

V/c c / 7-0+ ^A^o) J-oo / 

C k ,i(m,n;(3,a) = lim f— ^- /°° dv D m (P,v)* D n (a,v) h k>l {v)*) (179) 

\fc c / 7^0+ YA 1 ^oJ J-oo J 

where J = Ind [/i/^] +m + n and m + n > 2 is required in ( |1 77| ) . The resonance denominators 



in D m (P,v)* D n (a,v) are defined using our previous notation ( |111| ), ( |113| ), and ( |114| ); in 
particular, the pole locations are: 

atj = z + ijUj/k c j = 1, . . . , n (180) 

/3* = z* - i 7 IK j = l,...,m, (181) 

and the sets of poles are indicated with the notation a = (ai, . . . , a n ), (3 = (Pi, ... , p m ), 
v = (ui, . . . , u n ), and ( = . . . , ( m ). In using the notation above, e.g. for S(m, n; P, a), it 
is understood that the first m arguments after the semi-colon correspond to the poles in the 
lower half-plane D m (P,v)* and the remaining n arguments denote the poles in the upper 
half-plane D n (a,v). 

The latter two limits are related by 

C kjl (m, n; P, a) = e« [B k}l (n, m; a, P)}* (-l) m+n ^ (182) 

where the phase e 1 ^ is defined by 

e ieC*0 = i im ( ^M-) . (183) 
7 -o+ V A'i(zo) J 

Note that since A' k (zo) = d k (zo) for 7 > 0, this is the phase mentioned in the Introduction. 
From ( |171| ), the coefficients of the leading singularities are given by 
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bj = £ 2)i -i(0, 2; zq, zo) - S 0j -i(0, 2; z Q , z ) (184) 
+ S 2,i-i-2(0, 2; z , «o) - c h i B ,j-i-2{0, 2; z , zo)] 

1=0 

where is defined from ( |151| ) by the limit 

Cl , = lim ( 7 4{m) T M ) . (185) 

7^0+ v ' 



This expression is evaluated below, c.f. ( |2 10| ) . 

The properties of S(m, n; (3, a) in ( |177| ) are analyzed in Appendix B; the main conclusion 
needed here concerns the dependence on \x c . 

Lemma V.l For non-negative integers (m,n) such that m + n > 2, S(m,n; (3,a) has the 
form 

S{m, n; (3, a) = d{m, n; C, v) + {-l) m+n d{n, m; u, () e i?(M (186) 

where the real-valued functions d(m, n; (, v) are independent of Fq(v, yU c ). Thus S(m, n; a, (3) 
depends on fi c only through the phase exp(i£). 

Proof. See Appendix B. □ 

The analysis of Bkj(m, n; (3, a) follows the pattern for the calculation of hkj(v) in Table 
I. From the definition (|178| ), B fi(m, n; j3, a) and £> 2j o(^, n; (3, a) can be evaluated in terms 
of S(m, n; a, (3): 



B 0fi (m, n;p,a) = J2 S(m + 2,n + l; (3, f3 h z , a, z ) (187) 

n 

+ S(m + l,n + 2;(3,zo,a,ai,z ) 
i=i 

m 

S(m + 1, n + 2; (3, (3 U a', z ) (188) 

i=i 

n+l 

+ ^2 S{m^ n + 3; (3, a', a[, Zq) 



B 2 , (m,n;P,a) = - 



i=i 



i=l 



where a' = (a, z 2j0 ). This expression for B2,o( m , n ] (3, a) uses the identities: 

k 2 - 1 

\imk k {z Kl ) = —— (189) 

7->o+ k 
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from (|33|), and 



^ D* D 



r)(v,ii) = -6 m>0 S n0 . 



(190) 



-oo V — Z k ,l 

The latter identity is verified by noting that the integrand has index m + n — 1 soifm + n > 
we get zero, and when m = n = then the integral reduces to A 1 (^) — 1. From these 
expressions for £> ,o an d £> 2i0 their dependence on /i c is easily characterized. 

Lemma V.2 For all non-negative integers (m,n), Bo t0 (m, n; (3, a) and B2fi(m,n; j3,a) de- 
pend on fi c only through the phase exp(z£). 

Proof. 



This follows immediately from ( |187| ) - ( |188| ) and Lemma V.l. □ 



Recursion relations for the remaining limits Bkj(m,n; (3,a) and also for Ck,i in (|151| ) 
follow from the truncated relations in ( |172| ) - (|176|) ; this derivation is briefly summarized for 
k — 0,1, 2, and k > 3. These relations are somewhat complicated, but provide the basis for 
Theorem V.l below. 



1. B 0t i(m,n;(3,a) forl>l 



From the definition (|178|) and substituting for hoj(v) from ( |172| ) 



B r(m,n;P,a) 



m+n—1 



7 



lim . 

7 ^o+ \2 7 (/ + l)A / 1 (z ) 



dv D m (/3, v)* D n (a, v) I o l (v) 



191) 



Using ( |173|) the integral becomes 

dvD m (p,v)* D n (a,v) i i(v) 



(192) 



oo 

l-l 



■E(i+j')(ph+Ph; 

3=0 



dvD m ((3,v)* D n (a,v) h ,j(v) 



1-2 



h h-i - ^ c r^_! + Kj T i,i-j-2 

3=0 



CC 
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With the substitution 



d m D* D n D* D 

_ _ in* n \ — \ " m " _i_ \ " m n 



the limit in ( |191| ) can be rewritten as 



i=l 



v — a, 



(193) 



2(7 + l),Bo,z(m, a 

-E(!+i) 



(194) 



lim ^-^(pi-j+pt •_<) 



7~>0+ 
m+n m 



m+n— 1 



7 



4j+m+n+l ^oo 



+ (t) £. ln s 



+ 



m+n n 



dv ~ m 2 

oo V — Pi 



£ l im [1 TT ^— [ dv m n 



lim . 

7-0+ V Ai(zo) 



and then evaluated to obtain 

2(1 + l)B ,i(m, n; i 8 I a) = -5](l + j)(&f-j + 6?-,) S oj(m, n; /3, a) 

i=o 



(195) 



+ E ^ c i,2-i( m + !> n ! A Aj ") - + 1, ra; /?, A, a) 

i=i [ 

-c* ll _ l S{m + l,n+ l;/3,Pi, z ,a) + cy_i«S(m + 2,n;/3,^, z , a) 



Z-2 



+ E [ c i,i-i-2Cij(m + l,n;P,Pi,a) - c* 1> i_ j _ 2 B ltj (7n + l,n;(3,/3i 

3=0 
n ( 

+ 22 I Ci,i-i(m, n+l;(3, a, a { ) - E x ,i-\(m, n + l;(3,a, ai) 
i=i { 

-Cy_! S(m, n + 2; f3, z , a, a») + ci^-i <S(m + 1, n + 1; /3, z , a, "i) 

Z-2 

+ E [ci,/-i- 2 Cij(m,n + l;P,a,ai) - c* ll _ j _ 2 Bi t3 {m,n+ l;/9,a,Oi) 
i=o 



£ Bi j(m, n; /?, a) for I > 



From the definition (|178| 

B lt i(m,n; (3, a) 



m+n— 1 



lim 



7 



4Z+m+n+3 „oo 



a/ 7-0+ ^ Ai(zb) 

Using = h,i and ( p2|) for Fy, this can be rewritten as 



dv D m (/3, v)* D n (a, v) h lt i(v) . (196) 
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B h i(m, n; /?,«) = - 5(m, n + 1; /?, a, zy) 



(197) 



m+n 



lim 



7 



4i+m+n+3 



>A; c y 7^0+ y A'i(^o) ^-oo u - z 1)t j 

where we have used (|185|) . Using the expression for Ii t i(v) in (|174|) we obtain for the second 
term 



i 

k 
i-i 



m+n 



lim 



7' 



Al+m+n+'Z 



7-0+ \ A'^Zq) 



- D m (P,v)*D n (a,v)h ;l {vy 



-^[(2 + j)b l ^ j + (l+j)bt j 



j=0 



dv 

v - z u 

B 1;j (m,n + l;/3,a,zi,i) 



(19 



+ 



B 0> i(m, n + 2;(3, a, z 1)h z 1;l ) - B 2 ,i{m, n + 2;(3, a, z lh z 1)t ) 



i-i 



+ J2 ( c U - j -iB Qtj (m,n + 2;f3,a,z 1:l ,z 1 j) - c^^S^m, n + 2; (3, a, z u , z lit 

j=0 V 



-S(m,n + 2;P,a,z 1;h z ) 



B ,i{0, 2; z , z ) - B 2>1 (0, 2; z , z ) 



i-i 



3=0 



+ Ci,l-j-iBo,j{0, 2; z , Zq) - c^_ i _ 1 B 2 j(0, 2; z , Zo) 



+ ^2 \ B 0t i{m + l,n + 1; (3, pi,a, z h i) - B 2j j(m + 1, n + 1; /?, A, a, zy) 
i=i 



+ H ( c i,/-i-i' S o,j("T. + l,n + A>a>*i,i) - <?\,i-3-\®%i{ m + ™ + 1; A A, a > *U 



+ H I ^o,z( m ; + 2; /5, a, a h z X j) - B 2 ,i(m, n + 2;(3, a, a u z hl ) 
i=i 



i-i 



+ H ( c M-i-i^o,j(m,n + 2; /3, a, ai, - c* lj i_ j _ l B 2 ,j(m, n + 2;/3, a,^,^) 

3=0 V 



5. B 2> i(m, n; (3, a) for I > 1 



From the definition 



z \ 



m+n— 1 



lim 



7 



4/+m+ra+l 



B 2 ,i(m } n;f3 } a) = 7- , 

VA; C / 7-0+ \ Ai(«b) 

Using /i 2) i = R 2 {w 2 j) I 2) i, the integral becomes 



cfo D m (/3, v)* D n (a, v) h 2> i(v) . (199) 



dv D m (/3,v)* D n (a,v) h 2t i(v) 



—1 

2k r . 



- dv D* m D n I 2il (v) 

-oo v — Z 2 i 



(200) 



oc 



dv 



D* m D nV (v) r i 2>l (v' 



4A 2 (Z 2 ,0 ./-oo u - Z 2 ,Z J -00 v' - z 2;i 



dv' 



43 



and with (p9|) - (|T90|) we obtain 



1 ('Y^ (l\ 6m >° 6n A lim / 7 4 ' +m+n+1 f°° D* m D n h,i{v) 



(201) 



Let a' = (a,Z2 ) i) and define D n+ i(a' ,v) = D n (a,v)/(v — z 2) i), then using (|175|) the limit in 
( TO yields: 



m+n 



lim 



7 



4i+ra+n+l 



°° dv D m (p,vyD n {a,v)i 2 j(v) 



k c J - •()• \ A', (.:,,) ./-x 

-E[( 2 +j)6i-i+J^i 
i=o 



(202) 



B 2 ,j{m, n + l;/3,a') 

fll s 

+ Y, \ + 1, n + 1; A, «') - B 3) /_i(m + 1, n + 1; /?, A, a') 

i=i 

-ci jZ _i5(m + l,n + 2;(3, (3 h a', z ) 

i-2 x 

+ Y [ci,i-j-2Bij{m + l,n+ l;/3,/3i,a') - cl^^B^m + 1, n + 1; /3, A, a') I 

n+l , 

+ ^ j iB lji _ 1 (m, n + 2; /3, a', a-) - j6 3)/ _i(m, n + 2; /3, a', a-) 
i=i 

-ci j i_i«S(m, n + 3; /3, a', a-, z ) 

+ 12 [ci,i-j-2Bi,j(m, n + 2;(3, a', a-) - c^.^^- (m, n + 2; /3, a', a-) 



i=o 



^. Bk,i(m, n; (3, a) for k > 3, I > 



From the definition (|178|) 

m+n— 1 



B kt i(m,n;/3,a) = ( — 



lim 



7" 



7-0+ VAi(^o) 

where J = 2k + 4Z — 3 + m + n. Using fa^ = R k (w k> i) h,i 

D* m D n I k>l (v) 



dvD m (P,v)*D n (a,v)h k)l (v)\ (203) 



— i 



dv D m (/3,v)* D n (a,v)h k} i(v) = — 



dv 

oo V - Z k i 



dv 



(204) 



k 2 A k (z Kl ) J-oo V - Z Kl J-oo v' - z kji 



dv'' 



hi(v') 



and with (pgp - (|190D we obtain 
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B k ,i(m,n; (3, a) 



1 / i 



k \k. 



m+n 



j ^m,0 8n,o \ / 7 J f 00 dv Dm{l3,vY Dn{a,v)ikA V ) 



v - Z k i 



(205) 



Let a' = (a,Zk,i) and define D n+ x(a' ,v) = D n (a,v)/(v — z k ,i), then from (|176|) we find 



i \ m+n / 7 J 
*J 7™+ V A i(^o) 



Z-l 



3=0 



B kJ (m,n + 1; /?,<*') 
(206) 



/ft s 

+ { Bk-i,i(m + 1, n + 1; /?, A, «') - + 1, n + 1; /3, A, «') 

i=l 

l-l 1-2 . 

+ X! CM-i-i^fc-ijC 7 " + 1, n + 1; /3, A, a') - ^ ^,j-j-2^*+ij(»" + 1, n + 1; /3, A, a') \ 

3=0 3=0 ' 

+ { S fe-i,K m > n + 2 ; A «'> «D - 5fe +1) i_i(m, n + 2;f3, a', a-) j 

i=l J 

Z-l i-2 x 

+ Ci,i-j-iBk-i,j( m ' n + 2 i & a '> a - H c *,*-j-2#fc+ij( m > n + 2; /?, a', a-) I . 

« — n — n J 



5. Evaluation of Ck t i(fJL c ) 



From ( |151| ) and (|92"1), we have 



Ck,i(Vc) = lim 

7^0+ 



7 



(2fc+4Z-3+55 M ) 



(207) 



for > 1 and Z > 0. Since Ai(zk,i) = Ai(zo + i r ydk,i/k c ) = i'jdk,iA' 1 (zo)/k c + 0('y 2 ), as 7 — > + 
the prefactor gives 



~%k j k(- 



oh) 



k = 1 



\k 2 -l + A 1 (z k ,i)J 

where d\j = 2(1 + I). Thus for k — 1 
ci,i(/i c ) = - lim 



(208) 



-ik/k, 



k 2 - 



f + 0( 7 ) fc > 1 



7 



4Z+3 



7^0+ \d Xyl M x {Zq) J-oo V - ZlJ J ' 

and the right hand side is evaluated by setting m = n = in ( 198|) : 



(209) 
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(i){-§[ (2+j)6 - +(1+ *- 



^,,(0,1;^,) 



(210) 



+ 



B Ot i(0, 2; z lh ztj) - S 2 ,z(0, 2 ; z i,h Z\,i) 
i-i , 

+ Z ( c U-j-i#oj(0, 2; z^,) - cJ^.j-.jBajCO, 2; zj,,, z^)) 
i=o v 

B ,z(0, 2; 2b, «b) - S2,i(0, 2; z , Zo) 
i-i , 

+ Z ( c i,;-i-i^o,j(0, 2; z , 2b) - c^.^Ba^O, 2; z , ^ 
j=o v 



In deriving ( |210| ), the relation 5(0, 2; zij, z ) = 1 from ( |B10| ) in Appendix B has been used. 

The corresponding expression for c^i when k > 1 is much simpler. Setting m = n = in 
( TO and O, yields 



?A' 1 (^ c /fc c + ?0+) K? 



(211) 



from (207 ) for k > 2. Here u c /k c is the phase velocity at criticality; note also that when 
m = n = the arguments (j3, a) in Bk,i(m, n; f3, a) are irrelevant. 



C. Evaluation of the lowest order singularities b\ and 62 

It is instructive to calculate the cubic singularity b\ and the fifth order singularity 6 2 
from these integrated recursion relations. In the cubic case this merely recovers the previous 
result ( |168| ), but the calculation illustrates the formalism. The procedure is to apply the 
recursion relations until expressions involving S(m, n;/3,a) are obtained, then the results 
from Appendix B are used. 

1. Calculation of b\ 



For j — 1, the expression for bj in ( |184j ) yields 

bi (/i c ) = #2,0(0, 2; z , z ) - i3 ,o(0, 2; z , z ), 



(212) 



and from (|187 ) - ( 188|) one has 
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£>o,o(0, 2; zq, Zq) — 2 5(1, 4; zq, Zq, Zq, Zq, Zq) 



(213) 



^2,o(0, 2; zo, Zo) — «5(0, 5; Zo, Zq, 2:2,0, Zo, Zo) + 9^^' ^' Z2 ' ' ^ 2 ' ' Zo )' (214) 
In Appendix B the functions S{m, n; j3, a) are evaluated; applying these results yields 



^0,o(0, 2; Zq, Zq) — 2 5(1, 4; Zq, Zq, Zq, Zq, Zq) 



d{l, 4; 0, 0, 0, 0, 0) - d(A, 1; 0, 0, 0, 0, 0)e i? 

2 d(l, 4; 0,0, 0,0,0) 
1 

4' 



and it follows from (|B13| ) that £> 2j o(0, 2; z , Zq) = 0. Hence, as expected, b\ 



(215) 
■1/4 from 



2. Calculation of b 2 



The evaluation of 62 proceeds similarly but is considerably more laborious; I summarize 



the calculation below but omit the details. For j = 2, (|184|) yields 



b 2 {Hc) = #2,i(0, 2; zq, zq) - B 0jl (0, 2; z , z ) - Ci j0 B 0) o(0, 2; z , z ) 



(216) 



since B 2 ,a(0, 2; z , z ) = from (gig) and ( pT3|) . Setting / = in ( gI5D gives 



c i,o — 2 [^o,o(0, 2; z , z ) — B 0)0 (0, 2; z lj0 , z^q)] 



(217) 



With £>o,o(0, 2; z , Zo) = 1/4 from (|215|) and a similar evaluation yielding <6o,o(0, 2; zi j0 , z^o) 
1/32, this gives ci )0 = 7/4 3 and Q216| ) becomes 



b 2 (fi c ) = B 2) i(0, 2; z Q , z ) - #o,i(0, 2; z Q , z ) 



44; • 



(218) 



Consider i3 ,i in ( [21 8| ) first; from the /c = recursion relation 



B 0i1 (0,2;zo,Zq) = - 



-Bo,o(0, 2; z , z ) + 2 (V € 23i i0 (3, 0; z , z , z )* - Bi, o (0, 3; z , z , z ) 



+ 2 Ci^ <S(1, 3; Zq, Zq, Zq, Zq) 



(219) 
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and since c^o <S(1, 3; zq, zq, zq, zq) = (7/4 3 )(— 1/4), this simplifies to 



£o,i(0,2;z ,^o) 
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+ e J? 6i, (3, 0; z , z , z )* - B lfi (0, 3; z , z , z ) 



For 62,1 in (|218|) , from ( |201| ) and ( |202|) and using the results -62,0(0, 3; zq, zq, Z2,i) 
5(0, 5; P, a) = 0, one finds 

B 2 ,i(0,2;z ,z ) = - 



2 ( jBi i0 (0, 4; Zq, Zq, Z 2 ,l, Zq) ~ 63,0 (0, 4; Zq, Zq, Z 2 ,l, Zq) 
+6i,q(0, 4; Zq, Zq, Z 2 ,l, 22,1 ) _ 63,0(0, 4; Zq, Zq, Z 2 ,l, Z 2 \ 



(220) 
and 

(221) 



From the recursion for k = 3 in (|205|) , two terms in (|221|) vanish: 63,0(0, 4; z , z , z 2> i, Zq) = 
and -6 3 , (0, 4; z , z , z 2}1 , z 2) i) = 0, leaving 



6 2 ,i (0,2; 2:0,2:0) 



6i,o(0, 4; Zq, Zq, Z 2j \, Zq) + -6i,o(0, 4; Zq, Zq, Z 2) l, Z 2 .\) 



(222) 



The remaining terms in (|220|) and Q222| ) involving 61,0 are more tedious to evaluate from 
the k = 1 recursion relation in ( 1971) ; they have the following values 

15 e^ c) 



6i,o(3,0 ; Zq, Zq, Zq) 

6l,o(0,3 ; Zq, Zq, Zq) 
6l,o(0,4 ; Zq, Zq, Z 2) i, Zq) 

6l,o(0,4 J Zq, Zq, Z 2} i, Z 2j i) 



31 

2« 



2 6 



379 



(3 3 )(2 8 ) 
107 

WW) 



(223) 
(224) 
(225) 

(226) 



From these results ( 220 ) and (|222 ) give 

60,1(0,252:0,2:0) 



41 



9 



and from (|218 ) 



B 2!l (Q,2;zo,z ) = — 



(227) 
(228) 



(229) 



The most striking feature of this result is that, like the cubic coefficient b\, the fifth order 
coefficient does not depend on the critical equilibrium F (t>,/i c ). 
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D. Dependence on ji c : the role of 



The recursion relations in (PI , ( TO - QIH5D , f55p, (p7|) - (P5|) , ( M ) - (ffl), (|D|) - 
( j206|) , and (|210|) determine a closed set of equations for the coefficients {bj} and {ci /} and the 
functions {Bk,i} (also determined from (|182|) ). These relations are not independent 

of the critical equilibrium F (f,/i c ), but the dependence on the parameters /i c is entirely 
through a functional dependence on the phase expi£(/z c ). This implies that the higher order 
coefficients 63,64, . . . can only depend on F (v,fj, c ) through their dependence on expz£(/i c ); 
a more precise statement is given in Theorem V.l below. 

In order to describe the iterative procedure to be followed, it is helpful to organize this 
collection into nested subsets T>(N) as follows. For iV = let 



V(0) = {B 0fi ,B 2fi } 



(230) 



and for N > 1 define 



V(N) = 



{bj(Vc) I j 



N} 



{ciM \ l = 0,...,N-l} 



{Bo tl \l = 0,...,N} 



(231) 



{B ltl \l = 0,...,N-l} 

{B Kl I k = 2, . . . k max (N) and 1 = 0,..., (k max (N) - k)} 

where k max (N) = N + 2. In this notation it is understood that Bkj denotes the entire set 
of functions Bkj{m,n; (3, a) for all non-negative integers (m,n). Obviously T>(N — 1) is a 
subset of V(N) and as iV — > 00 all coefficients and functions mentioned above are included 
in V(N). 
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Lemma V.3 For N = 1,2, . . the recursion relations determine all elements ofT>(N) in 
terms of the elements ofV(N — 1), the functions {S(m, n; j3, a)}, and the phase e^. 



Proof. 



This is readily verified by inspection of the recursion relations, and tak- 
ing into account the identity Ci s i(m, n; (3, a) = (— l) m + n_1 Bi s i(n, m; a, (3)* expi£ 
from flTgp . □ 



The recursion relations lead to our main result concerning the dependence of the 7 — > + 
limit on the underlying critical equilibrium F Q (v,[i c ). 

Theorem V.l For < N < 00, the elements of T>(N) depend on the critical parameters 
[i c only through a functional dependence on the phase . In particular, for 1 < j < 00 
there exist functions Qj (z) , satisfying 

Q j (z)* = Q j {z% (232) 

such that 

b j { fJlc )=Q j (e i ^). (233) 

Each function Qj is universal in the sense that it is independent of k c and Fq(v,{i c ). The 
other elements of V(N) can also be similarly expressed as functions of the phase exp(i£) 
which satisfy fl232|) . 

Proof. 

The proof is by induction. 

1. Lemma V.l shows that S(m,n; /3,a) depends on /i c only through a func- 
tional dependence on and the functions involved are simple polynomials 
with real coefficients. Specifically S(m, n; (3, a) = s(e^) where 

s(z) = d(m, n- 0, a) + (-l) m+n d(n, m; a, (3) z. (234) 
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It then follows immediately from ( |187| ) - (|187|) for i3 ,o and i3 2 ,o that these 
functions can also be expressed as polynomials of e^. Since d(m, n\ /3, a) is 
real- valued, these polynomials satisfy ( [23 2| ). This proves the theorem for 
P(0). 

2. For the induction step assume that the theorem holds for T>(N — 1) for 
some fixed N > 1, and consider the recursion relations which determine 
T>(N) from T>(N — 1) and {S(m, n; /3, a)}. Each of these relations is a sum 
of terms and most of these terms fit the following description: the term 
consists of either a single element of T>(N — 1), a product of two elements 
of D(N — 1), or a product of an element of T>(N — 1) with S(m, n; (3, a) (for 
some (m, n)). In each case the term is multiplied by a real coefficient which 
is independent of \i c and k c . The exceptions to this description are discussed 
below. For each term covered by this description if the individual elements 
of T>(N — 1) appearing in the term are functions of fi c only through e^, then 
the entire term has this property. In addition, if the individual functions 
satisfy ( |232| ), then their product will also. Thus all terms, describable in 
this way, will depend on /x c only through a functional dependence on 
and this function will satisfy (|232|) . Hence these terms preserve for T>(N) 
the functional dependence on /i c assumed for V(N — 1). 

3. The recursion relations also contain two types of terms that differ from the 
above description. First, there are terms which are the complex conjugate 
of an element of T>(N — 1), or which are the product of an element of 
T>(N — 1) with the complex conjugate of a second element of T>(N — 1). In 
either case the coefficient of the term is real and independent of \i c and k c . 
For example in ( |184| ) one finds the term c\ l £> 2 j_z_2(0, 2; zq, Zq). Since by 
assumption there are functions fi(z) and f2{z) such that c\ t i = /i(e^) and 
£>2,j-z-2(0, 2; zq, z ) = /^(e^), this term leads to a functional dependence 
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cj,i B 2 j-i- 2 (0,2]Zo,zo) = /(e*) where f(z) = h{l/z) f 2 (z). Therefore this 
type of term also preserves for for T>(N) the functional dependence on \i c 
assumed for T>(N — 1). 

4. The second type of exceptional term involves Cij and only arises in the 
recursion relation for B j in ( |195| ). Typical examples are C\j,-x(m + 



1, n\ (3, /3j, a) and Ci^_j_ 2 Ci,j( m + 1, n; j3, /3 i7 a). Using the identity (|18 
these terms can be re-expressed in terms of Bij as: 

Ci,,_i(m + 1, n; (3, ft, a) = (-l) m+n e « B 1)M (m + 1, n; (3, A, «)* 
ci^-aCuCm + l,^/?,^,^) = (-l) n ^e«B u (m+l,n; i 9,ft I a)*. 

In each case the essential difference from terms already discussed is that the 
coefficient, (— l) m+n e^, is complex. However since this complex coefficient 
is obviously a function of the phase as well, it does not change the conclu- 
sion: these terms also preserve for for T>(N) the functional dependence on 
fi c assumed for T>(N — 1). Hence if the theorem holds for T>(N — 1), then 
it holds T>(N). Since the conclusion has been verified above for T>(0), by 
induction the theorem holds for all N. 



□ 



When the equilibrium has reflection symmetry Fo(v,fi) = Fo(—v,fi) then p(a,/j,) is real 
and hence bj = Qj(e 1 ^) must be real also. One can check that in the case of reflection 
symmetry A[(z ) is pure imaginary and expi£ = —1. Thus the reality condition (|232|) 
ensures that in this case bj = Qj(—1) is real as expected. 

Note that the explicit calculation of b\ and b 2 in the previous Section reveals that the 



first two functions in ( 233 ) are constant: 



Q 1 (z) = ~ (235) 
13 

Q2(z) = - A , (236) 
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thus non-trivial dependence on expz£ can only arise in bj for j > 3. It is natural to wonder 
if the Qj are simply constants to all orders. Although this cannot be categorically ruled 
out, an inspection of the recursion relations does not seem encouraging. The evaluation of 
( |236| ) involves a crucial and seemingly accidental cancellation that eliminates all the terms 
depending on expz£. If such a cancellation persists at still higher order, then detecting it 
will require a deeper understanding of the integrated recursion relations. 



E. Implications of Theorem V.l 



Theorem V.l implies that as 7 — > + the dynamics represented by our amplitude expan- 
sions can depend on the critical equilibrium F (v,fi c ) only through a functional dependence 
on derivative of the dielectric function evaluated at the phase velocity of the critical linear 
mode u c /k c , 



A[(uj c /k c )= lim Ai(z ) 
7^0+ 



p.v. r dv ^ v ^\ 

00 (v-uj c /k c ) 



9ti 

i7i—(u c /k c ,fi c ). 



(237) 



This may also be written in terms of the phase (|183|), 



A[(Lu c /k c ) = |Ai(w c /A; c )|exp(-i£/2); 



(238) 



in fact for many features it is only the phase exp(i£) that matters and not the magnitude. 
For example, the asymptotic form of the amplitude equation ( |16C| ) 

dr 



dr 



.2/ 



(239) 



depends on F (t>,/i c ) through bj, and from (|233|) this can be re-expressed as 

dr 



dr 



l + ^fReQ^e^)]^' 



(240) 



where Qj{z) itself is a universal function, i.e. independent of Fq. For the electric field in 
( |167|) , the Fourier component at the wavelength of the unstable mode is 



\k c E l {t)\= 1 2 r{r) l+Y\r{r 



(241) 
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where 

oo 

r; = E c i^ 2/ - (242) 

1=0 

Since according to Theorem V.l, the coefficients {cij} also depend on F (v, /i c ) only through 
a universal functional dependence on exp(z£), this component of E(x, t) is predicted to have 
an asymptotic dynamics that is determined by exp(i£). For the other wavelengths, one finds 
that there is an overall factor of A'^Uc/kc) in Ej,(t) (cf. (|2 1 1|) ) , so that \k c Ek(t) / A[(uj c / k c )\ 
is determined by exp(i£). 



VI. THE DISTRIBUTION FUNCTION: ASYMPTOTIC BEHAVIOR 

The evolving distribution function F(x, v, t) = F (v, p) + f u (x, v, t), rewritten using 
([5l"l) and (|53"D , takes the form 



F(x,v,t) = F (v,fi) + ah (v,a) (243) 



+ 



pit) (Mv) + <rhi(v, a)) e^" 9 ^ + ^ p k h k (v, a )e tk ^ x - 9 ^ + 



cc 

k=2 



where the phase and amplitude variables p(t)e %e( ^ = A(t) have been used and a = \A\ 2 . 
For fixed 7 > 0, our analysis of the mode amplitude dynamics leads to the equations 



- Q1571 ), but the long time evolution of p(t) is difficult to predict since the higher order 
nonlinear terms are not negligible; this difficulty persists even when 7 — ► + as shown in 
( |160|) . However if p(t) asymptotically approaches a constant value as t — > 00, 

pit) -> Poo, (244) 



then 6(t) must settle down to a fixed frequency d9/dt ujoo according to ( [L57| ). In this 



event, the form of F(x : v : t) in (|243j) must approach a travelling wave moving at a constant 



wave velocity u^/fcc. Hence the limiting behavior in fl244j ) implies the flow on the unstable 
manifold asymptotically approaches a travelling wave such as a Bernstein-Greene-Kruskal 
mode, pi This asymptotic state is necessarily periodic in time due to the periodic boundary 



conditions. 
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The asymptotic form of F(x, v , t) as 7 — > + can be analyzed from (|243| ) without needing 
to know the long time evolution. In terms of the rescaled amplitude p(t) = 7 2 r(7t), ( |243| ) 
becomes 



F(x,v,t) = F (v,fi) + r 7 h (v,a) 



(245) 



+ 



r { 1 2 i, c {v) + 1 & r 2 h l {v,a)) e^ ^ + £ 7 2fc r%(^, a)e^-^» + 



cc 



fc=2 



For /raec? v ^ v„, as 7 — > + , this implies 



F(x,v,t) - F (v,n) 



r 



r ij c (v)e i{kcX - m) + ccl + C(7 2 ); 



(246) 



thus away from the phase velocity the correction at the wavelength of the unstable mode 



is dominant and is given by the critical eigenf unction. At v — v p , ( 245 ) must be examined 
more closely since the functions hk(v, a) are singular when 7 = 0. 



A. Asymptotic behavior near v = v p 



Our main interest is in the structure of F(x, v, t) near v = v p as 7 — » + . It is necessary 
to extract the singular behavior of hk(v p ,a) and balance it against the explicit factors of 
7 in ( 245 ). The key idea, pointed out by Larsen [41], is to appropriately magnify the 
neighborhood of ^(7) using a rescaled velocity coordinate for 7 > 



u 



7 



'-(.v -v p (i)). 



(247) 



The motivation for this definition is the resulting factorization of a resonant denominator: 



1 



v - z k l 



7 



(248) 



_u - i(l + d k ,i) 

in which the 7 _1 singularity has been extracted and the remaining function of u is nonsingular 
as 7 — > + . 

Consider the effect of this coordinate change on the linear eigenfunction (|3"4]). Note that 
since d v F (u c /k c , /i c ) = 0, 



55 



( 1 \ dF 

T](v p + nfu/ k c ,n) = - ( — j -Q^Ml) + W M(t)) 



(249) 



" 7 Uf. 



<9 2 F n 



<9t> 2 



dv p (0) u 



<9 2 F 0/ ,d/i(0) <n . . 



__7 



dvdfi 



11 

- + Re[A' 1 (uj c /k c )} + 0( 1 ) 



g?7 



(250) 



where in the last step the first order solutions for ^(7) and /i( 7 ) from (|A13|) in Appendix 
A have been used. This motivates the definitions 



rj(u,fjt)= \^Jr)(v p + ju/k c ,iJ,), 
if) c (u) = i>c{v p + iu/k c ) 



(251) 
(252) 



so that 



-r)(u,n) 



u — l 



(253) 



Thus the eigenfunction has a non-singular limit as 7 — > 0; from ( 250 ) - fl251|) 

[Re[A[(uj c /k c )]+u/k 2 c } 



lim ^c(u) 
7-^0+ 



lu — i) 



(254) 



The generalization of these definitions for the nonlinear theory is easily done for any 
function characterized by an index as defined in Section IV. Let Q(v,fi) have index Ind [Q] 
then define Q(u,fj,) 



<?K/i)^ 7 1+Ind Wgivr + Tu/kc,!*)] 



(255) 



this is consistent with ( |252| ) since the eigenfunction has index —1. It is not hard to see that 
Q(u,fi) will be non-singular as 7 — > + . In particular this definition accomplishes the goal 
of extracting the singularities in hki' 



^(«) = 7 1+Ind [hk ' l] h k:l (v p + 1U /k c ). 



(256) 



Applying this to hk(v, a) gives 
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hk{v p + 7«/ k c , a) = ^ + 7«/ &c) o"' 



Z=0 

oo 



1 



l+Ind 



2/ 



E7 4i M h k Au)r 
1=0 \ ' / 

2fc-2+4(<5 M +5 fc ,i) oo 
Z=0 



i 



(257) 



so if J = 2k — 2 + 4(Sk,o + Sk t i) then r ) J h k {vp + ^u/k cl a) is non-singular as 7 — > + . Formalize 
this observation in the definition 



h k (u,r 2 ) = (7) 



2fc-2+4(5 fc , +5 M ) 



h k (v p + 'ju/k c ,'y 4 r 2 ) = h k>l (u)r 



,21 



(258) 



1=0 



The left hand side of ( |246| ) can now be evaluated sX v = v p + ^u/k c in terms of i/j c (u 
and hk(u, r 2 ), 

[F(x, v p + ju/ k c , t) - F (v p + ju/ k c , //)] 



where 



g(x, u, t, li) = r 2 h (u, r 2 ) + 



(259) 



r (ip c (u) + t hi(u, r )) e 



2\\ A{k x-e{t)) 



+ r ' k h k (u,r 2 )e 



2\ ik(k c x-9(t)) 



+ CC 



k=2 



(260) 



As 7 — > + , this yields a nonsingular expression g(x, u, t, /i c ) for the distribution function in 
the neighborhood of the phase velocity. In contrast to (|246|) , here all wavelengths contribute 
to the leading correction at O^ 2 ). 

It is natural to consider whether g(x,u,t, fi c ) is in some sense universal, or equivalently 
to ask how does g(x,u,t, li c ) depend on F (v,fj, c )l Obviously there is a trivial dependence 
on the linear frequency u through the factors expik(k c x — 9(t)) which can be suppressed by 
considering the Fourier coefficients: 



\g (u,t,Li c )\ = r 2 h (u,r 2 ) 
\gi(u,t,fi c )\ =r 4> c (u) +r 2 h 1 (u,r 2 ) 



(261) 
(262) 
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and for k > 2 



\g k (u,t,n c )\ = r k h k (u,r 2 ) 



(263) 



these depend on F (v,fi c ) through ip c (u) and the functions hk(u,r 2 ). At \i = fi c , the eigen- 
function ( |254p depends on the critical equilibrium only through the derivative of the dielec- 
tric function A[(u! c /k c ), and the dependence of hk(u,r 2 ) on F (f,/i c ) can be investigated by 
analyzing the series coefficients hk,i{u) in (|258|) . The explicit forms for /io,o(w) and h,2fl{u) 
follow from (H - (H) 



^2,0 (u) = - 



d 




T)(U, 




9 




du 


|> 


"*) 


(« + *)_ 


du 


_{u + i)_ 




H 








(u 


— i 




6(u — 





J 



(264) 
(265) 



as 7 — > + , these expressions depend on \i c only through if) c (u) and therefore depend on 
F (v,fi c ) only through A[(u) c /k c ). 

The remaining coefficients are determined iteratively with recursion relations that follow 
by making the change of variable v = v p + 7u/fc c in the truncated relations in ( |172|) - ( |1 76|) 
and allowing for the asymptotic behaviors in ( |151| ) and ( |152| ). Let 

Ind[ ^ ] [4, i K+7«/y = 7 IlldK ' 



hA u ) = 7 
then for k = 



h,i{v p + Tu/k c ) +0(7) 



(266) 



hoAu) 



Iq,i{ u ) 
2(1 + 



(267) 



with 



J o,i( u ) 



ES(l+j)(6i-i + 6?- 4 -)M«) 



/ = 



z > 



(268) 



&l,J-l(u) - C U-l^c(M) + E5=0 ^lj(w)ci,i-j-2 



cc 
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For k > 1, the general relation 

hk,i(u) = 7 1+Ind {hkA [(Rk(w k>l ) I k>l ) (v p + 1U /k c )\ ; 



(269) 



applies, although the evaluation of the resolvent depends on whether k — 1 or k > 1. For 
k = 1 



7 



i+Ind [ ftM ] [( i2l ( Ti;i|| )/ li| )( t;p + 7U / fcc )] 



i[u - i(l + tZi,i)] 



+ 0(7) 
(270) 



where 



i-i 



h,iW = "E [(2 + j)bi-j + (1+M-j] hiA" 



(271) 



j=0 

9 



i-i 



+i du i ~ + ^ (V( M ) c ^-;-i - h 2,j( u ) c h-j-i) 



3=0 



iip c {u) f 00 du' 



i-i 



A[(uj c /k c ) J~oo (u'-i) 2 

+0(7), 

and for k > 2 

y+Ind [hk,i\ [{Rk{wKl) {Vp + 7m/A;c) ] 



MO - M U ') + E (/lOj(«') C l,i-J-l - ^j( U ') C l,i-j-l) 
3=0 



Ik.lM 



ik[u - i(l + d K i)} 



+ 0(7) 



where 



*2,j( U ) = ' 



- E l -J [(2 + j)bi-j + jbQ h 2tj (u) 

+i lk {Ki-i( u ) - h,i-i(u) + ^ c («)ci,i-i 

+ £5=0 - M^K*-^] } + 0(7) 



Z = 



I > 



and 



h,i( u ) = "EP + + 6 *-j] frfcjH 

9 



«-2 



(272) 



(273) 



(274) 



+2 <j /ifc-i,i(w) - /ifc+i,i-i(«) + E h k-i,j{u)ci,i-j-i - h k+i,j( u ) c *i,i-j-2 
+0(7) 
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for k > 3. 

Inspection of these recursion relations leads to a simple characterization of the \i c depen- 
dence of g(x, u, t, fi c ). 



Theorem VI. 1 The Fourier components \gk(u,t } /i c )| of g(x,u } t } fi c ) in ( |261| ) - ( |263| ) de- 
pend on the critical equilibrium F (v, fi c ) only through a functional dependence on A' x {oj c /k^) 
where 



M x {u c /k c ) = lirn Ai(zb). 

7— >0+ 



(275) 



Proof. 



It suffices to examine the dependence of i[> c (u) and hk(u,r 2 ) on Fo(v,[i c ). The 
/i c dependence of ip c (u) is through Re [h! x {u c /k^\ from (|254|) , and the de- 

(J,=fl c 

pendence of hk(u,r 2 ) = hk t i{u)r 21 will be inferred from the coefficients hk,i{u). 



1. By inspection from (|264|) - ([26 5|) , the lowest order coefficients in the orga- 
nization of Table I also depend on F (v,fi c ) through Re [A[(u) c /k c )]: 



ho )0 (u) 
h 2 ,o(u) 



d_ 

du 



Re[Ai(o; c /A; c )] +u/k, 



d 



u 2 + l 



Re[Ai(o; c /fc c )] + 



u — i 



2{u — i) du 

2. The remaining coefficients {hk,i(u)} are generated from {ip c , ho,o, ^2,0} using 
the recursion relations ( |267[ ) - ( |274| ) which depend explicitly on fi c through 
A'^ujc/kc), fj(u,fi c ), and the coefficients {bj} and {c X) i}. By Theorem V.l, 
the coefficients depend on /i c only through the phase 

= A[(u c /k c y 

" KMK) ' 

and from ( |250| ) - ( |251j ) fj(u, fi c ) depends on fi c through Re [A[(u c /k c )]. Thus 
by induction, the property that {^ c , h Ot0 , ^2,0} depend on \i c only through 
A[(u) c /k c ) is passed on to all the coefficients {hk,i( u )}- 



(276) 
(277) 



(278) 



□ 
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VII. DISCUSSION 



The approach followed in this paper is patterned on the well established methods of 
center manifold reduction and normal form analysis which have proved quite powerful in 
analyzing bifurcations in dissipative systems. However for the Vlasov equation the methods 
must be adapted. For a dissipative system undergoing a Hopf bifurcation, one can reduce the 
problem to a finite-dimensional submanifold for all /x sufficiently close to /x c , and for /x > /x c 
this corresponds to the finite-dimensional unstable manifold associated with the equilibrium. 
By contrast, in the Vlasov instability there is no analogous possiblility of reducing to a finite- 
dimensional submanifold in a full neighborhood of /x = fi c since the critical eigenvalues first 
appear at /x = /x c embedded in the continuum and then emerge for /x > /x c . It is only on one 
side of criticality (/x > /x c ) that the linear spectrum indicates finite-dimensional submanifolds 
should exist and offer a useful parallel to the dissipative case. 

The center manifold and normal form analysis for Hopf bifurcation can be rephrased in 
terms of the unstable manifold by deriving the dynamics on the unstable manifold and then 
taking the 7 — > + limit of the resulting vector field to obtain normal form equations for the 
amplitude A(t). As pointed out in the Introduction, the scaling A(t) = ypy ' r(^t)e~ iel ^ then 
yields asymptotic equations [ f42f 

— = rR H (r,ri (279) 

GST 

HO 

- = 00 + 0( 7 2 ) (280) 

where 

R H (r, /x) = 1 + [ ai (ji c ) + (9( 7 )]r 2 + 0( 7 2 ). (281) 

This result reveals the correct scaling ^7 for the nonlinear oscillation and shows that the 
expansion Q281| ) for the amplitude dynamics can be truncated after the lowest order nonlinear 
term. After truncation, the remaining dependence on the underlying critical equilibrium is 
expressed entirely through the lowest order nonlinear coefficient ai(/x c ) which is a calculable 
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function of /i c . This truncation also allows the time-asymptotic behavior to be determined 
since the solutions of Rh = are then easily found. 

Our application of this same procedure to the instability of an electrostatic wave in a 
collsionless plasma leads to equations like ( |279| ) - ( |280| ) with several qualitative differences. 
First, the scaling of the mode amplitude required to obtain nonsingular asymptotic equations 
is quite different: A{t) = 7 2 r^^e" 10 ^' where now 

dv 

— = rR{r,fi) (282) 

GST 
rlf) 

- = u + 0( 7 ) (283) 



with 



R(r, 11) = 1 + Y, [R e Qj(e ia " c) ) r 2j + 0(7)]. (284) 



Secondly, the expansion for the amplitude equation ( 284 ) cannot be truncated. This means 
that the solutions to R(r, /x c ) = are not readily found and consequently ( |282| ) makes no 
obvious prediction for the long time behavior. However despite the greater complexity of 
R(r,fjL c ), the dependence on \i c can be characterized as arising only through a functional 
dependence on the phase e 1 ^^ which is in turn determined from the derivative of the 
dielectric function. 

This characterization of R(r, fi c ) (and similar statements for the electric field and dis- 
tribution function) makes testable predictions about instabilities driven by physically very 
different distributions. If e l ^ c > is fixed, then any variations in densities or temperatures 
characterizing Fo(v,[i c ) do not affect the evolution of r(r). For example, a beam-plasma 
instability (complex A) and a two-stream instability (real A), compared at a common value 
of e 1 ^, have identical amplitude equations up to 0(7) corrections. Examination of this pre- 
diction through numerical solution of the Vlasov equation will be undertaken in a future 
paper. 

A final important difference with the dissipative case concerns the generality of the 
analysis. In standard center manifold reduction, the local attractivity of the submanifold 

62 



ensures that the time-asymptotic behavior of any initial condition near the equilibrium can 
be reliably predicted from the evolution on the submanifold. For Vlasov, or more generally 
in a Hamiltonian bifurcation, when there are neutral modes in addition to the critical modes, 
then there may be little or no correlation between the evolution of an initial condition on the 
unstable manifold and the evolution of an arbitrary initial condition. Nevertheless numerical 
studies of the one mode instability observe the trapping scaling in the saturation amplitude 
of the electric field as a quite robust phenomenon; there is no evidence that the initial 
condition must be carefully chosen. This indicates that some features of the evolution on 
the unstable manifold have a wider validity, but it is not clear what selects these features. 
One appealing conjecture is that for initial conditions near F , but not on the unstable 
manifold, the electric field E(x,t) evolves asymptotically towards the electric field E u (x,t) 
associated with a solution on the unstable manifold. This could occur without requiring a 
corresponding asymptotic behavior in F(x, v, t) and would suffice to explain the robustness 
of the trapping scaling. This picture seems difficult to investigate analytically but can be 
tested through numerical experiments. 

A basic feature of the dynamics expected from an autonomous one-dimensional flow such 
as dr/dr = R(r, /i c ) is the absence of oscillations, provided R is differentiable in r (or at 
least Lipschitz continuous). This absence of oscillations implies that if r(r) approaches an 
asymptotic limit r(r) — ► r^as r — > oo, then the approach will be monotonic from below. 
This is indeed found in Hopf bifurcation where -R#(r, aO is simply a quadratic polynomial in 
r. However in numerical studies of the one mode instability such monotonic relaxation in the 
time-asymptotic regime is not observed, rather one finds the familiar trapping oscillations 
in the electric field. 

The explanation of trapping oscillations within the setting of unstable manifold dynamics 
may be related to the survival of an infinite sum of terms in ( |284| ) as 7 — > + . The fact that 



higher order terms in r are not higher order in 7 is a marked contrast with the much simpler 
limit found in the Hopf normal form (|281[ ), and is directly related to the fact that the critical 
eigenvalues merge with the continuous spectrum as 7 — > + . In exactly solvable models, 
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where neutral modes also introduce singularities into the amplitude equation, I have shown 
that the unstable manifold can develop a spiral structure which persists as 7 — > + and this 



spiral allows the flow to approach an asymptotic limit through a decaying oscillation. |43 
This spiral structure is illustrated in Fig. 4. 

When such a spiral is present, then describing the dynamics on the manifold via a 
mapping H from the unstable subspace (|45l) yields a vector field on E u with branch point 
singularities at the points where the flow moves from one branch of the spiral to the next. 
In the solvable examples, the fact that R(r, /x c ) has a branch point within the domain of 
flow implies that the higher order terms in the expansion of R(r, /x c ) remain essential to 
the dynamics even as 7 — ■> + . Supposing that the amplitude dynamics for the one mode 
problem has a similar structure, then R(r, fi c ) would have a branch point at r&, the turning 
point of the spiral. As the mode grows, the increase of r(r) to r& would signal the onset of 
trapping oscillations with the passage of the trajectory to the next branch of the unstable 
manifold. 

Note that a trajectory will reach such a branch point node in finite time, unlike the more 
familiar situation of a node where the vector field is differentiable and the approach time 
is infinite. In addition, the loss of smoothness at r = rj, introduces the lack of uniqueness 
needed by the solution to pass through the branch point. Finally it should be clear that 
such a spiral structure would present a significant obstacle to using the power series (|284|) 
to determine the time- asymptotic amplitude r^. 

Instabilities in other systems, including ideal shear flows [jl~4-4~6||, solitary waves |47|-p9|, 
bubble clouds f50fl , and globally-coupled populations of oscillators also exhibit key 

features of this problem, most notably that the unstable modes correspond to eigenvalues 
emerging from a neutral continuum at onset. In the case of ideal shear flows similar singu- 
larities arise in the amplitude equations for the unstable modes (7 -3 at cubic order) p6| ; 
by contrast, in globally coupled phase models for the onset of synchronized behavior in a 
population of oscillators the critical eigenvalues emerge from the continuum at the onset of 
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instability but the amplitude equations are nonsingular and ^7 scaling is found (at least 
in the best understood case of a real eigenvalue). |37|j51| , [53| This difference in the nonlinear 
behavior seems noteworthy since the linear dynamics of the oscillator model is qualitatively 



similar to Vlasov although apparently lacking a Hamiltonian structure. |52| In the models of 
solitary waves and bubble clouds, the scaling behavior in the weakly unstable regime has not 
been investigated. It would be interesting to determine if singularities arise in the amplitude 
equations for the unstable modes in these problems with corresponding implications for the 
scaling behavior of the nonlinear states. 

The study of these novel bifurcations from a unified viewpoint is just beginning, and it is 
not yet clear how to abstract the essential features required to produce singular expansions 
and unusual scaling behavior. 
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APPENDIX A: PARAMETRIZED ROOTS OF THE DIELECTRIC FUNCTION 

For 7 > 0, the roots of the dielectric function determine the eigenvalues of the unstable 
modes. From ( |29"1) this requires the phase velocity v p (j) and parameter ^(7) to satisfy 
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A 1 (^ o ( 7 ),/i( 7 )) = l--/( 7 ) = 



(Al) 



where 



°° dvd v F (v,n(~f)) 



3 v - v 9 (i) - ij/k c 

For simplicity, /i is taken to be a single parameter. 

From ( |Al| ) these functions can be calculated perturbatively for small 7; let 



Kl) =/i c + /ii7 + C(7 2 ), 



(A2) 



(A3) 
(A4) 



and also expand /( 7 ): 



J( 7 ) = /(0)+7^(0) + O( 7 2 ) 



Then ( |A1| ) implies 



1(0) = k 



dl 

g?7 



(0)=0 



(A5) 



(A6) 
(A7) 



and so forth; the coefficients (i>i,//i) are determined by (|A7|) . Define v = v p (j) + u/k c and 
take u as the variable of integration in /, 

du 8Fq 



then 



£«» = p.v. r - 

CX7 J-oo u 



00 u — vy dv 



(i \ d 2 F c) z F 

k~ + V ~dv^ ^ c// kc + M// kc ' ^ + ^ dv¥~ ^ c// kc + M// k " ^ 



(u p (7) + w,/z( 7 )), 



(A8) 



<9 2 F n 



dvdfi 



k r ' V <9t> 2 



d 2 F 0/ „ _ «9 2 F n 



■ (u c / k c , He) + ^1 - - (a> c / fe c , /i c ) 
<7t> a/i 



(A9) 



Since the real and imaginary parts must vanish separately in ( |A7[ ), there are two equations 
for Vi and H\. 
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1 



oo du d 2 F 
P.V. / —-(u} c /k c + u/k c ,n c ) 

oo U OV z 



P.V. 



00 du 

-00 u 



<9 2 .F <9 2 .F 
vi -q^- {oj c / k c + uj k c , /i c ) + Vig^Q^ (uc/ k c + u/ k c , /i c ) 



7T 



d 2 F [ 
dv 2 



u c /k c ,fi c ) 



(A10) 



(All) 



With 



A[(uj c /k c ) = - 



kl 



P.V. 



00 du d 2 F 



uj k c + uj k c , /i c ) 



-00 ' u dv 2 

from (|237|) , these equations may be rewritten more simply as 



d 2 F 1 

+ i-K-^-(u c /k C7 ii c )\ (A12) 



P.V. 





r^F . 




[ dv 2 {U 


POO 


du d 2 F 


1 — 00 


u 


dv 2 




P.V. 



dvdfi 
uj c /k c + u/k c ,fi c ) 



Vc/k c ,Hc 



k c Re[A[(uj c /k c )] (A13) 



-oo-^ ( " cAc + MAc ' /ic) 



-ym[Al(w c /fc c )]. (A14) 



APPENDIX B: EVALUATION OF SINGULAR INTEGRALS 



Let D = 1 and for n > define 

«) = t ry ^ r (Bl) 

(v - ati)(v -a 2 )---{v- a n ) 

where a = (ai, . . . , a n ). For non-negative integers (m, n) such that m + n > 2, we consider 



limits of the following type (c.f. (|177|) ): 

m+n—2 



S(m,n;P,a) = i^—j Inn f ^ dv D m (j3, v)* D n (a, v) T](v, //) J , (B2) 



7^0+ y A[(zo) J -00 
where the poles of the integrand 



atj = z + i^fj/kc j = 1, . . . , n (B3) 

/5; = ^-*70 /k c j = l,...,m (B4) 
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lie along the vertical line Re v = v p . The non-negative constants Vj > and Q > are 
assumed to be independent of F for all j; in particular they are independent of 7. The limit 
S(m, n\ ft a) does depend on F as determined below. 

Several general relations can be noted immediately; first interchanging the order of the 
arguments gives the simple identity 

S(n, m; a, (5) = (-l) m+n e ^ c) S{m, n; ft a)* (B5) 

where 

= Hm I (B6) 
is the phase defined previously in ( |183| ). Thus it is sufficient to evaluate (|B2|) for m < n. 



Secondly, by expanding the denominator with partial fractions, the limit S(m, n;j3,a) can 
be expressed in terms of the limits for m + n — 1 . Ifm>l then 

[S{m- l,n;P",a) -<S(m- l,n;(3',a)) 
S (m, n; /?, a) = — (B7) 

Sm—l Sm 

can be used, if n > 1 then 

[<S(m,n- l;(3,a") - S(m,n - 1; (3,a')) 

S(m, n;p, a) = (B8) 

v n — v n-\ 

applies, and if m > 1 and n > 1 then 

[«S(m- l,n;^',a) -S(m,n- l;P,a')] 
S(m, n; (3, a) = — — (B9) 

can be used. In these recursion relations, the primed arguments are defined by a' = 
(q;i, . . . , a n _i), «" = («!,..., a n _ 2 , a n ), /5' = (ft, . . . , /3 TO _i), and = (ft, . . . , /3 m _ 2 , /3 m ). 



1. Proof of Lemma V.l 

These general properties allow S(m, n\ (3, a) to be calculated recursively from results for 
m + n = 2, 3. For m + n = 2 the integral is straightforward to evaluate and the limits are 
given by 

68 



«S(0, 2; ax, (X2) = 1 



S(2,0; 

5(l,l;A,ai) 



z/x + Ci e* 



2 + vx + Ci 

For m + n > 2, when either m = or n = 0, the integral is nonsingular as 7 
results in (|B10| ) and (|B11| ) generalize easily 



<S(0,n; a) = <5„ j2 . 
S(m,0;/3) = 5 m , 2 e 4 «. 

A similar evaluation for m + n = 3 gives two more limits 

(2 + Ci)-Ge* 



(BIO) 
(Bll) 
(B12) 

+ , so the 

(B13) 
(B14) 



5(l,2;/5i,ai,a 2 ) 
5(2, 1; ft, &,<*!) 



(2 + z/ 1 + Ci)(2 + z/ 2 + Ci) 
1/1 - (2 + i/i)e* 



(2 + z/ 1 + Ci)(2 + z/ 1 + C 2 )' 
By inspection of ( |B10| ) - ( |B16| ), the results for m + n = 2, 3 have the following form 



(B15) 
(B16) 



S(m, n; ft a) = d(m, n; C, + (-l) m+n d(n, m; z/, C) e i? 



where £ = (Ci, • • • , Cm), v — (vx, ■ ■ ■ , v n ), and the functions d(m, n; (, v) are given by 



(B17) 



d{0,2;ux,iy 2 ) = 1 
d(2,0;Ci,C2) = 

d(i,i;Ci,^) = 2 + ^ +Ci 



(B18) 
(B19) 
(B20) 



for m + n = 2 and 



d(l,2;(x,vx, U2) 
d(2,l;Ci,C2^i) 



2 + 



(2 + z/ 1 + Ci)(2 + z/ 2 + Ci) 
^1 



(2 + z/ 1 + Ci)(2 + z/ 1 + C 2 ) 
for m + n = 3. 

The representation in (|B17p in fact holds for all m + n > 2. From ( B13 ) and 



(B21) 
(B22) 
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d(0,n;v) =6 n , 2 (B23) 
d{m, 0; C) = 6 m>2 . (B24) 



More generally, when m and n are both non-zero, inserting ( |B17] ) into ( |B7| ) - ( |B~9| ) yields 
the corresponding recursion relations for d(m, n; (, v): 

d(m, n; (, u) = ^ - ^ C>) " ^ - ~ " " > (B25) 





- C, v) - 


• d(m, n — 






2 + v n 






d(m, 


n-i;C/)- 


- d(m, n — 














-l,n;C'»- 


- d{m — 1 


n; C, f) 



d(m, n; C, ^) = 2 L22 -^ — (B26) 



d(m, n; C, f) = t 37 (B27) 

Srrt— 1 W 

where (' = (Ci, ■ ■ ■ , Cm-i), = (^1, ■ ■ • , ^n-i), C" = (G, • • • , Cm-2, Cm), and 1/' = 
(1/1, . . . , z/ n _2, u n ). In ( |B26| ) and ( |B27|) there is no singularity when the denominator vanishes 
because the numerator also vanishes; more convenient forms for these recursions are given 
below. 

Thus the recursive evaluation of S(m, n; /3, a) reduces to the recursive evaluation of 
d(m, n;C>^)- Note that the functions d(m, n; C, v) are real-valued and universal in the 
sense of being independent of the equilibrium F (v,fi c ). The limit S(m, n; f3, a) depends 
on Fq(v,/j, c ) only through the phase e^. This completes the proof of Lemma V.l. 

2. Evaluation of 62 

In the remainder of this Appendix, some results for the functions d(m, n; (, z/ ) that are 
useful in the calculation of 62 are briefly summarized. The case when m = or n = is 
trivial, and from ( |B21| ) and the recursion relation in ( |B25| ), it is easy to show that for n > 2 



m£?Z?v (B28) 

Thus the functions are explicitly known unless both m and n are greater than 1 so that 
m + n > 4. 
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Henceforth assume that both m and n are non-zero and that m + n > 3. In (P2 1|) - 
( |B22| ), note that d(2, 1; Ci, C2, ^1) = 2; 1/1 — 2, 2 + Ci, 2 + £ 2 ); more generally such a relation 
also holds for m + n > 3: 
d(m, n; Ci, ■ ■ ■ , Cm, v u . . . , u n ) = (-l) m+n_1 d(7i, m; v x - 2, . . . , v n - 2, 2 + &, . . . , 2 + Q- 

(B29) 

This identity is readily verified by induction from the recursion relations for d(m, n; C, v) in 
( |B25|) - (|B27| ). Hence it is sufficient to calculate d(m,n; C, v) for m < n. 

From the definitions it is clear that d(m, n; (, v) must be symmetric under interchange 
of the m arguments (d, . . . , Cm) and also under interchange of the n arguments (u±, . . . , u n ). 
It is convenient to have a notation that makes this more explicit and which is manifestly 
nonsingular when some of these arguments coincide and the denominators in ( |B26| ) and 
( |B27| ) vanish. For m > 1 and n > 1, define N(m, n\ (, v) by the formula 

N(m,n;(, v) 



d(m,n;(,v) 



Y\U nr=i(2 + n + 



(B30) 



so that 



JV(l,l;Ci,J>i) = fi 
iV(l,2;Ci,^i,^) = 2 + Ci 
iV(2,l;Ci,C2^i) = v x . 



(B31) 
(B32) 
(B33) 



From ( |B25|) the corresponding recursion relation for N(m, n; C, v) when m and n are greater 
than 1 and m + n > 4 is 



JV(m, n; C, 1/) = iV(m - 1, n; f ) 



'n-l 



n( 2 +^+Cr> 



i=l 



iV(m,n-l;C,0 



m— 1 



H(2 + v n + (,•) 
(B34) 



and from (B29|) we have the identity 
N(m, n; &, . . . , Cm, vi, ■ ■ ■ , u n ) = (-l) m+n - x iV(n, m; ^ - 2, . . . , v n - 2, 2 + Ci, • • • , 2 + Cm) 

(B35) 
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when m + n > 3. 

Clearly the functions N(m, n; (, v) are also symmetric under interchange of the m argu- 
ments (Ci, • • • , Cm) and also under interchange of the n arguments (y\, . . . , v n ) but this is not 
manifest in ( |B34j ) . Some functions N(m,n;£,i/), useful in the evaluation of b 2 , are listed 
below showing explicitly the interchange symmetry; in these formulas the notation P n cor- 
responds to the symmetric polynomials defined by Pi(x, y) = x + y, Pz(x, y) = x 2 + xy + y 2 , 
P 3 (x, y) = x 3 + x 2 y + xy 2 + y 3 , and Pi(x, y) = x 4 + x 3 y + x 2 y 2 + xy 3 + y 4 . 

N{2, 2; C, v) = -v x v 2 + (2 + Ci)(2 + C2) (B36) 
N{2, 3; C, v) = v x v 2 v z - (2 + Ci)(2 + C2) {Pi(2 + Ci, 2 + C2) + vi + "2 + ^3} (B37) 
N{2, 4; C, u) = -v x v 2 v z v A + (2 + Ci)(2 + C2) {p 2 (2 + Ci, 2 + C2) + Pi (2 + Ci, 2 + C2) [fi + ^2 + ^3 + u 4 

+v\(v2 + v-i + va) + ^2(^3 + ^4) + ^3^4 1 (B38) 



N(2, 5; C, v) = ^ 2 *W 5 ~ (2 + Ci)(2 + C2) [p 3 (2 + Ci, 2 + C2) 

+P 2 (2 + Ci,2 + C 2 ) 



(B39) 



Vl + V2 + "3 + V<i + ^5 

+Px(2 + Ci, 2 + C2) [^1(^2 + ^3 + ^4 + v 5 ) + v 2 (y z + + u b ) 

+ Z/ 3 (l/ 4 + ^5) + ^4^5 



+V1 



^2(^3 + ^4 + V5) + ^3(^4 + ^5) + ^4^5) 



+ ^2 



1/3(1/4 + z/ 5 ) + z/ 4 i/ 5 



+ ^3^5 



and 



-(2 + Ci)(2 + C 2 )(2 + Cs) {(2 + Ci)(2 + C2 + 2 + Cs) + (2 + Ca)(2 + Cs) (B40) 

+ (6 + Ci + C2 + CsX^i + V2 + V3) 
+v\ + v\ + v\ + v\{yi + z/ 3 ) + v 2 vz I 

+^1^3 |^l(^2 + ^3) + ^3 + (6 + Cl + C2 + C3)(^l + V2 + ^3) + (2 + Cl) 2 

+ (2 + C2) 2 + (2 + C3) 2 + (2 + Ci)(2 + C2 + 2 + Cs) + (2 + C 2 )(2 + Cs)} 
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N(3, 4; C, u) = (2 + Ci)(2 + C 2 )(2 + Cs) {(2 + Ci) 2 [(2 + C2) 2 + (2 + Cs) 2 ] (B41) 

+ (2 + C 2 ) 2 (2 + Cs) 2 + (2 + Ci)(2 + C 2 )(2 + Cs)[6 + Ci + C2 + Cs] 

+(1/1 + v 2 + 1/3 + ^4) [(2 + Ci) 2 (4 + C2 + Cs) + (2 + C 2 ) 2 (4 + Ci + Cs) 

+ (2 + Cs) 2 (4 + Ci + C2) + 2(2 + Ci)(2 + C 2 )(2 + Cs) 
+ (z/ 2 + „l + u j + ^ [ (2 + Ci)(2 + G + 2 + Cg) + (2 + C2)(2 + Cs) " 

+[1/1(1/2 + ^3 + ^4) + ^2(^3 + ^4) + ^3^4] (6 + Ci + C2 + C3) 2 

+ (6 + Ci + C2 + C3) [^1(^2 + vz + ^4) + ^2(^1 + ^3 + ^4) + ^1(^1 + V2 + ^4) 

+v\{y\ + l/ 2 + I/3) + 2(i/ii/ 2 l/3 + l/ll/ 2 l/4 + ^1^3^4 + 1^1/3^4) 
+l/l(l/2 + 1/3 + 1/4) + l/|(l/3 + 1/4) + ^3^4 + I/? (^2^3 + ^4 + ^3^4) 
+l/f (1/1I/3 + I/1I/4 + ^3^4) + ^3(^1 ^ 2 + ^1^4 + ^2^4) 
+I/|(l/ll/ 2 + I/1I/3 + l/ 2 l/ 3 ) + 3z/iz/ 2 z/ 3 z/ 4 | 

-1/HW4 {(2 + Ci) 3 + (2 + C2) 3 + (2 + Cs) 3 + (2 + Ci) 2 (4 + C2 + Cs) 

+ (2 + C 2 ) 2 (4 + Ci + Cs) + (2 + Cs) 2 (4 + Ci + C 2 ) + 2(2 + Ci)(2 + C 2 )(2 + Cs) 
+(1/1 + i/ 2 + 1/3 + 1/4) [(2 + Ci) 2 + (2 + C2) 2 + (2 + Cs) 2 

+ (2 + Ci)(4 + C2 + Cs) + (2 + C 2 )(2 + Cs)' 
+ (6 + Ci + C2 + C3) [^1(^2 + ^3 + ^4) + ^2(^3 + ^4) + ^4 

+ ^l(^2^3 + ^2^4 + ^3^4) + ^2^3^4| 
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FIGURES 

FIG. 1. Linear stability boundary for the beam-plasma instability described by a Lorentzian 
plasma and a Lorentzian beam as in ( |37|) with L = 2ir, n = 0.8, A = 0.3, and u p = 0.0. The 
instability occurs at the longest wavelength corresponding to k = 1; modes with k = 2,3, .. . are 
always stable. The linear spectra for (a) stable equilibria, (b) the critical equilibrium and (c) 
unstable equilibria are illustrated in Fig. 2. 

FIG. 2. Spectrum of C near criticality for the beam-plasma instability of Fig. 1. (a) The 
subcritical spectrum contains only the continuous spectrum and a (degenerate) eigenvalue at zero 
which is related to the degenerate Hamiltonian structure of the Vlasov equation. |34|] The continuous 
spectrum coincides with the imaginary axis but is slightly thickened for ease of visualization, (b) 
At criticality, the conjugate pair of eigenvalues, (A, A*) with A = 7 — iu, appears for the first 
time embedded in the continuous spectrum, (c) The supercritical spectrum shows a quadruplet of 
eigenvalues in addition to the continuum and zero eigenvalue. 

FIG. 3. Local geometry of the unstable manifold; the equilibrium Fq is at the origin. 

FIG. 4. Conjectured spiral structure in the global unstable manifold shown in cross section with 
the 6 coordinate suppressed. The turning points of the spiral correspond to branch points, e.g. 
r = r;,, in the mapping functions describing the manifold. The time-asymptotic state at r = is 
not on the branch of the manifold connected to the equilibrium at r = 0. 
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TABLES 

TABLE I. Order of calculation of ht,i(v) and pj from ij) c (v). The flow of calculation of the 
hkj(v) is indicated by moving downward. From ip c (v), /io,o an d /i2,o can be calculated and then 
pi determined; h± t o and Z^o are calculated next from {p±, /io,o, ^2,0} and then /io,i and /i2,i can be 
evaluated. This then determines P2, and so forth. For N > 2, pjy requires prior calculation of hk,i 
forO</c<iV+landO</<iV-/c + l- 2(5 kfi + <5 M ). 





k = k = 1 k = 2 k = 3 fc = 4 fc = 5 k = 6 


Po 




Pi 


fy),0 - ^2,0 


P2 


^1,0 ^3,0 

h h 
"0,1 "2,1 


P3 


h\ : l /l4,0 
^0,2 ^2,2 ^3,1 


P4 


^1,2 ^5,0 
^0,3 ^2,3 ^3,2 ft-4,1 


P5 


^1,3 ^6,0 
^0,4 ^2,4 ^3,3 ^4,2 ^5,1 
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